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'  An  Investigation  of  the  interaction  of  electrostatic  waves  launched 
by  a  slow  wave  structure  with  a  magnetized  plasma  Is  made.  The  charac¬ 
teristics  of  the  electrostatic  waves  and  the  electron  dynamics  are  studied 
experimentally. 

Of  primary  experimental  interest  is  the  measurement  of  the  electron 
energy  distribution  and  the  rf-induced  electron  flux  along  the  background 
magnetic  field.  This  interest  is  motivated  by  a  need  for  a  more  complete 
understanding  of  interaction  of  plasma  with  a  slow  wave  electrostatic 
field  which  is  of  importance  for  rf-heating  and  rf  dc  current  drive  in 
plasmas .  .• - 

Electrostatic  waves  are  launched  from  a  slow  wave  structure  near  the 
low  density  periphery  of  a  magnetized  plasma  column  that  possesses  density 
gradients  in  both  the  radial  and  axial  directions.  The  axial  inhomogen¬ 
eity  is  found  to  significantly  complicate  the  picture  of  wave  propagation. 

It  is  found  the  waves  launched  along  the  background  magnetic  field  into 
a  region  of  Increasing  electron  density  can  experience  a  reflection  as 
the  wave  propagates  into  the  plasma  interior.  Significant  increases  in 
the  effective  electron  temperature  and  rf-induced  axial  electron  flux 
near  the  plasma  surface  are  observed  using  an  electrostatic  energy  analy¬ 
zer  and  current  probes.  The  lower-hybrid  resonance  layer  as  determined 
by  the  observations  of  large  radial  phase  changes  using  rf  probes  in 
this  region  correspond  approximately  with  the  radial  location  of  Increased 
electron  temperature  and  rf-induced  electron  flux. 

The  increase  in  effective  electron  temperature  can  be  attributed 
to  the  coherent  oscillations  of  the  resonant  electrons  in  the  wave 
field  as  derived  from  the  quasillnear  theory.  The  rf-induced  electron 
flux  can  be  interpreted  as  the  nonlinear  ponderomotive  force  on  electrons 
due  to  the  parallel  gradient  of  the  electric  field  energy.  The  experi¬ 
mental  results  are  in  qualitative  agreement  with  theory. 
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Chapter  1 
INTRODUCTION 


The  interaction  of  rf  waves  with  plasmas  has  been  extensively 
studied  especially  in  recent  times  because  of  their  applicability  to 
rf  current  drive  and  plasma  heating.  Techniques  involving  the  coup¬ 
ling  of  energy  to  lower  hybrid  waves  have  received  considerable  atten¬ 
tion  in  realizing  these  goals.  In  thermonuclear  devices  the  lower  hy¬ 
brid  frequency  is  in  the  low  gigahertz  range  where  the  high  power  re¬ 
quired  for  heating  and  current  drive  is  readily  available  using  con¬ 
ventional  high  power  microwave  sources. 

Wong^  has  observed  current  generation  by  unidirectional  electron 

2  3 

plasma  waves  and  McWilliams,  et.al.  *  have  studied  rf  driven  current 
by  collisionless  and  collisional  damping  of  lower  hybrid  waves.  In 
addition,  current  drive  using  lower  hybrid  waves  has  been  observed  by 
Bernabei,  et.al.  in  the  PLT  tokamak.  General  studies  of  rf  current 
drive  have  been  made  by  Bhadra,  et.al. ^  and  Fisch  and  Karney^.  Exten¬ 
sive  studies  have  been  made  in  the  literature  of  the  excitation  and 

7-13 

properties  of  lower  hybrid  waves  in  plasmas 

In  this  report  we  consider  an  experimental  study  of  the  in¬ 
teraction  of  plasma  with  electrostatic  waves  launched  by  an  external 

slow  wave  structure.  Electrostatic  waves  are  launched  along  the  back- 

2  “ilc  2 

ground  magnetic  field  with  either  an  e  2  or  e  2  dependence  at  the 
low  density  periphery  of  the  plasma  column.  The  plasma  possesses  density 
gradients  in  both  the  radial  and  axial  direction  which  complicates  the 
description  of  wave  propagation  in  the  plasma.  Experimental  investiga¬ 
tions  of  both  particle  and  wave  dynamics  under  the  influence  of  these 
externally  launched  slow  electrostatic  waves  are  examined. 

A  detailed  description  of  the  experimental  setup  is  given  in  chapter 
2.  The  plasma  production  scheme  and  the  system  used  to  excite  slow  elec¬ 
trostatic  waves  is  discussed  as  well  as  the  supporting  equipment  required 
for  the  experiment.  Chapter  3  describes  a  variety  of  diagnostics  used 
in  the  experiments  for  measurements  of  the  properties  of  the  plasma  as 
well  as  for  the  rf  waves.  Where  applicable,  the  basic  theories  of  the 
diagnostic  systems  are  discussed  so  that  a  reliable  interpretation  of 


experimental  data  can  be  made.  Detailed  experimental  results  are 
presented  in  chapter  4.  First,  a  description  of  the  background  plasma 
is  made  and  secondly  the  results  are  presented  from  the  measurements 
when  the  plasma  was  excited  by  the  slow  electrostatic  wave.  Measure¬ 
ments  include;  radial  and  axial  profiles  of  electron  density  and  temper¬ 
ature,  radial  and  azimuthal  profiles  of  axial  electron  current,  density, 
axial  electron  energy  distribution  function  as  determined  from  an  elec¬ 
trostatic  energy  analyzer,  and  amplitude  and  phase  measurements  of  the 
rf  potential  in  the  plasma.  In  chapter  5  a  brief  discussion  of  the  wave 
structure  in  a  plasma  with  both  axial  and  transverse  inhomogeneities  is 
made  using  the  cold,  magnetized  plasma  dielectric  tensor.  Also  a  dis¬ 
cussion  is  made  of  the  interaction  of  the  electrostatic  waves  with  the 
plasma  electrons.  An  interpretation  of  experimental  results  is  al? '  dis¬ 
cussed  in  this  chapter  with  suggestions  for  future  work. 
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.Chapter  2 
EXPERIMENTAL  SETUP 

To  understand  the  dynamics  of  a  plasma,  it  is  important  to  be 
familiar  with  the  method  used  to  produce  the  plasma  as  well  as  the 
diagnostic  techniques  used  to  obtain  experimental  data.  In  this 
chapter,  the  method  and  associated  systems  used  for  plasma  produc¬ 
tion  will  be  discussed.  Also,  the  system  used  to  excite  low  phase 
velocity  (u/k^<<c)  electrostatic  waves  for  the  rf  experiments  will 
be  discussed.  Chapter  3  will  be  devoted  to  the  various  plasma  diag¬ 
nostic  and  electronic  systems  used  in  the  experiments. 

2.1  General  Description 

The  plasma  is  produced  and  maintained  by  electron  cyclotron  reso¬ 
nance  heating  (ECRH)  at  2.45  GHz  in  a  magnetic  mirror  field  at  one  end 
of  a  linear  device.  The  plasma  drifts  out  of  this  coupler  region  into 
a  drift  region  with  a  uniform  axial  magnetic  field.  Measurements  are 
made  in  the  relatively  quiescent  drift  region.  A  general  pictorial 
representation  of  the  device  is  shown  in  Pig.  2-1.  Fig.  2-1  does  not 
show  the  slow  wave  rf  coupling  structure.  The  slow  wave  structure  and 
the  associated  rf  system  will  be  discussed  in  section  2.6. 

-4 

The  discharge  is  typically  operated  at  a  pressure  of  2  x  10  Torr 

with  Argon  as  the  filling  gas.  The  plasma  is  weakly  ionized  (typical- 

10  -3 

ly  0.5%)  with  electron  densities  on  the  order  of  10  cm  at  the  cen¬ 
ter  of  the  plasma  column.  The  electron  temperature  is  generally  about 
4  eV. 

2.2  Vacuum  System 

The  vacuum  system  is  composed  of  a  6-inch  high  speed  diffusion 
pump  backed  by  a  rotary  mechanical  pump.  A  Varian  model  336,  6-inch 
water  cooled  baffle  is  mounted  between  the  diffusion  pump  and  the  slid¬ 
ing  valve  to  prevent  backstr earning  of  pump  contaminants  into  the  vacuum 

chamber.  The  diffusion  pump  is  a  National  Research  Corporation  model 

-3  -9 

VHS6  with  an  operating  range  of  10  Torr  to  10  Torr  and  a  pumping 
speed  greater  than  2000  liters/sec.  The  mechanical  backing  pump  i 
Duo  Seal  model  1397  rotary  pump. 


Ionization  Gauge  -  J  rcas  Feed 


Basic  Experimental  Setup 
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The  mechanical  forepump  maintains  a  system  pressure  of  approxi- 
-3 

mately  10  Torr.  To  prepare  the  machine  for  plasma  production,  the 
diffusion  pump  is  activated  and  evacuates  the  system  to  10  ^  Torr  base 
pressure.  Argon  gas  is  then  leaked  into  the  chamber  using  a  flow  gauge 
to  obtain  an  equilibrium  flow  at  2  *  10  Torr  operating  pressure. 

-3 

Pressures  below  10  Torr  are  measured  by  a  Varian  model  563  Bay- 
ard-Alpert  ionization  gauge,  controlled  by  a  Varian  model  843  ratio- 
matic  ionization  gauge  control.  The  Varian  843  unit  also  contains  two 

_3 

thermocouple  gauge  monitors.  Pressures  above  10  Torr  are  measured  by 
a  Varian  model  531  thermocouple  gauge  tube. 

The  vacuum  chamber  which  basically  consists  of  three  main  sections, 
is  15.24  cm  in  diameter  and  has  a  total  length  of  approximately  1.2 
meters.  The  first  section  houses  the  microwave  energy  coupler  and  is 
constructed  of  aluminum  to  provide  adequate  shielding  of  the  microwave 
coupler  as  well  as  acting  as  the  reference  electrode  for  the  plasma  ex¬ 
periments.  Its  large  structure  which  is  externally  water  cooled  also 
dissipates  the  large  quantities  of  heat  generated  in  the  microwave  en¬ 
ergy  coupler.  It  is  approximately  11.9  cm  in  length.  The  second  sec¬ 
tion  of  the  vacuum  chamber  is  a  pyrex  4-port  tee  section  (see  Pig.  2-1). 
This  section  is  45.7  cm  long.  The  third  section  is  a  61  cm  long 
straight  pyrex  section.  This  last  section  has  fourteen  1.27  cm  diam¬ 
eter  ports  at  various  locations  for  plasma  diagnostics  and  the  power 
feeds  for  the  slow  wave  structure.  Pig.  2.2  shows  the  locations  of 
these  ports. 

2.3  Magnet  System 

The  axial  magnetic  field  is  provided  by  two  sets  of  Princeton  L-2 
water  cooled  coils.  One  set  of  ten  coils  establishes  the  magnetic 
mirror  field  used  in  the  ECRH  coupler  region.  These  ten  coils  are 
series  connected  to  a  Sorensen  model  DCR40-125  dc  power  supply  capable 
of  delivering  up  to  150A  at  a  maximum  power  of  7.5  kW.  The  power  supply  is 
operated  at  96A  establishing  the  ECRH  condition  for  plasma  production. 
The  second  set  of  coils  establishes  the  uniform  axial  magnetic  field  in 
the  drift  region.  These  ten  coils  are  series  connected  to  a  HuMac  dc 
power  supply  capable  of  delivering  uo  to  600A  at  a  maximum  cower  of  180 


kW.  The  maximum  allowable  current  through  the  L-2  coil  is  300A.  A 
uniform  magnetic  field  of  approximately  1.8  kG  results  if  the  drift 
coils  are  operated  at  the  300A  limit. 

Fig.  2-2  shows  the  location  of  the  L-2  coils  in  reference  to  the 
vacuum  chamber.  Fig.  2-3  shows  the  axial  magnetic  mirror  field  pro¬ 
file  per  unit  mirror  field  current  as  function  of  axial  position.  Also 
shown  in  Fig.  2-3  is  the  axial  drift  magnetic  field  profile  per  unit 
drift  field  current  as  a  function  of  axial  position.  Fig.  2-4  shows 
the  total  axial  magnetic  field  profile  for  a  mirror  field  current  of 
96A  and  a  drift  field  current  of  150A.  Most  of  the  experiments  were 
performed  using  these  magnetic  field  parameters. 

The  plasma  column  acquired  maximum  radius  at  the  region  of  mini¬ 
mum  magnetic  field,  i.e.,  near  the  region  of  the  high  vacuum  port. 

The  magnetic  field  is  uniform  over  the  circular  cross  section  of  the 
plasma  column. 

2.4  Microwave  Power  System 

The  microwave  system  used  for  ECRH  is  shown  in  Fig.  2-5.  The 
microwave  power  is  produced  by  an  Amperex  model  YJ1160  magnetron  cap¬ 
able  of  delivering  up  to  2.5  kW  cw  at  2.45  GHz  + 25  MHz.  The  anode  of 
the  magnetron  is  water  cooled  and  the  cathode  radiator  is  air  cooled. 
The  output  power  of  the  magnetron  is  proportional  to  the  beam  current 
as  shown  in  Fig.  2-6.  The  magnetron  operates  at  a  cathode  voltage  of 
approximately  -4.6  kV.  To  maintain  a  stable  microwave  power  output,  a 
high  voltage  beam  current  regulator  is  used  in  conjunction  with  the 
high  voltage  dc  power  supply  to  deliver  a  well  regulated  beam  current 
to  the  magnetron. 

The  output  coupling  of  the  magnetron  is  a  1  5/8"  rigid  coaxial 
line  with  an  Impedance  of  50fl.  The  magnetron  output  is  connected  to 
a  tapered  coaxial  line  which  reduces  the  1  5/8"  coaxial  line  to  a  7/8" 
coaxial  line  while  maintaining  a  line  impedance  of  50G.  The  7/8"  co¬ 
axial  line  is  coupled  to  a  WR-284  S-band  waveguide  through  a  coax-wave¬ 
guide  adapter.  The  microwave  energy  is  then  fed  through  an  isolator 
used  to  protect 'the  magnetron  against  excessive  reflected  power.  The 
isolator  can  dissipate  reflected  microwave  sower  up  to  750  W  average 
power.  The  output  of  the  isolator  is  connected  to  a  dual  directional 


Fig.  2-2  Vacuum  Chamber,  Field  Coil  Location,  Diagnostic 
Port  Locations  and  RF  Slow  Wave  Structure  Coupling  Ports. 


Fig.  2-5  ECRH  Microwave  System 


Fig.  2-6  Magnetron  Output  Characteristic. 
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coupler.  The  low  power  branches  on  the  directional  coupler  are  used  to 
monitor  incident  and  reflected  microwave  power  as  well  as  to  detect  ex¬ 
cessive  incident  or  reflected  microwave  power.  The  detection  circuits 
shut  down  the  high  voltage  power  supply  if  the  incident  or  reflected 
power  exceeds  750  W.  Incident  and  reflected  microwave  powers  are  mon¬ 
itored  by  an  HP  model  431B  power  meter  measured  through  59.9  dB  and 
60.6  dB  of  attenuation  respectively.  An  HP  model  536A  frequency  meter 
is  used  to  determine  the  microwave  frequency.  A  coaxial  switch  is  used 
to  select  which  directional  coupler  branch  is  connected  to  the  power 
meter.  The  main  branch  of  the  directional  coupler  is  connected  to  a 
7/8"  coaxial  line  through  a  waveguide-coax  adapter.  The  7/8"  coaxial 
line  is  brought  into  the  vacuum  chamber  through  a  7/8"  coaxial  pressure 
seal  and  tapered  to  the  microwave  energy  coupler  used  to  generate  the 
plasma.  The  microwave  energy  coupler  will  be  discussed  in  the  next  sec¬ 
tion.  A  thin  teflon  waveguide  section  1.59  mm  thick  is  inserted  be¬ 
tween  the  last  two  waveguide  sections  to  isolate  the  reference  electrode 
from  ground  (see  Fig.  2-1). 

2.5  Microwave  Energy  Coupler 

ECRH  is  achieved  by  coupling  electromagnetic  energy  into  the  plasma 
at  the  electron  cyclotron  frequency  using  a  Lisitano  The 

Lisitano  coil  has  been  used  previously  in  this  laboratory  by  Chrisner  ^ 
in  a  study  of  electron  cyclotron  resonance  plasmas  and  by  Kudyan  in 
a  study  of  transverse  instabilities  in  a  hollow  plasma  column.  The 
microwave  energy  coupler  is  shown  in  Fig.  2-7  and  the  physical  dimen¬ 
sions  of  the  coupler  are  shown  in  Fig.  2-8.  The  overall  assembly  of 
the  microwave  energy  coupler  is  shown  in  Fig.  2-9.  As  shown  in  Figs. 

2-7  and  2-8  the  energy  coupler  consists  of  an  open  ended  chin  copper 
cylinder  with  16  axially  oriented  slots  uniformly  distributed  over  the 
cylinder.  These  axial  slots  are  alternately  connected  together  by  short 
azimuthal  slots  at  each  end  of  the  cylinder  forming  one  continuous  slot 
around  the  cylinder.  The  first  slot  is  connected  by  a  fin-line  adapter 
to  the  coaxial  line  through  a  boron-nitride  vacuum  seal.  Boron  nitride 
was  chosen  for  the  vacuum  seal  because  of  its  ability  to  withstand  the 


:ig.  2-9  Overall  Assembly  of  Microwave  Energy  Coupler 
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large  quantities  of  heat  generated  by  the  discharge.  The  last  slot 
is  terminated  in  a  short  circuit.  The  length  of  each  axial  slot  is 
a  half  wavelength  long  at  Che  applied  microwave  frequency  (2.45  GHz). 
Without  a  plasma,  the  microwave  energy  is  guided  along  the  slots  to 
Che  shorted  termination  where  it  is  reflected  establishing  an  azi¬ 
muthal  electric  field  within  the  cylinder  with  a  maximum  standing 
wave  voltage  in  the  central  cross  section  of  the  cylinder.  The  elec¬ 
tric  field  configuration  is  shown  in  Fig.  2-10. 

After  a  plasma  is  established,  microwave  energy  is  absorbed  by 
ECRH  as  the  wave  propagates  along  the  slots  and  the  reflected 
power  is  determined  by  the  state  of  the  plasma.  In  turn,  the  state 
of  the  plasma  is  determined  by  the  coupler  region  magnetic  field, 
the  drift  magnetic  field,  and  the  neutral  gas  pressure. 

2.6  RF  System  and  Slow  Wave  Structure 

Fig,  2-11  is  a  schematic  diagram  of  the  3.25  MHz  rf  system  jjsed 
to  excite  low  phase  velocity  electrostatic  waves  in  the  plasma.  RF 
power  at  3.25  MHz  is  produced  by  an  rf  transmitter  capable  of  deliver¬ 
ing  up  to  5  kW  cw.  The  output  power  amplifier  consists  of  two  parallel 
connected,  grid  driven  Eimac  750  TL  power  triodes  operated  class  C. 

The  power  amplifier  output  is  connected  to  a  pi-L  network  which  tunes 
the  transmitter  freauency  and  adjusts  the  Dower  transfer  from  the  trans¬ 
mitter  to  the  load  impedance.  The  transmitter  output  is  connected  to 
a  5  kW,  wideband  (2MHz  -  30  MHz),  1  :  1,  ferrite  core  balanced  trans¬ 
former  which  is  in  turn  connected  to  two  50fJ  (RG-8/u)  delav  lines 
which  phases  adjacent  elements  of  the  slow  wave  structure  90°  apart 
as  shown  in  Figs.  2-11  and  2-12.  Two  power  meters  on  the  output  of 
the  transmitter  monitor  the  output  power  of  the  transmitter  and  the 
reflected  power  from  the  rf  coupling  system.  The  other  ends  of  the 
two  slow  wave  structure  delay  lines  are  terminated  by  50G  loads  each 
capable  of  dissipating  2.5  kW.  These  50Q  terminations  eliminate  re¬ 
flections  on  the  slow  wave  structure.  The  rf  load  currents  are  moni¬ 
tored  on  an  oscilloscope  or  an  rf  spectrum  analyzer  by  two  model  411 


Scope  or 


Fig.  2-11  Schematic  of  Low  Frequency  RF  System  and  Slow  Wave  Structure 


Fig.  2-12  Pictorial  Representation  of  Slow  Wave  Structure 
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Pearson  wideband  (1  Hz  -  33  MHz)  current  monitors  with  a  50Q  output 
impedance.  The  current  monitors  produce  an  output  voltage  of  100  mV 
into  an  open  circuit  with  an  inout  current  of  1A.  The  rf  spectrum 
analyzer  consists  of  an  HP  141S  oscilloscope  with  two  spectrum  analy¬ 
zer  plug-in  units;  an  HP  8553L  RF  unit  and  an  HP  8552A  IF  unit.  The 
spectrum  analyzer  is  capable  of  displaying  rf  spectra  from  1  kHz  to 
110  MHz.  Fig.  2-13  shows  the  rf  spectrum  of  the  transmitter  output. 
The  rf  output  and  terminating  loads  on  the  delay  lines  can  be  inter- 

X  1c  Z  •>  t  V  2 

changed  so  either  an  e  2  or  an  e  z  excitation  can  be  established 
on  the  slow  wave  structure. 

The  copper  coupling  plates  have  a  radius  of  5.1  cm,  a  length  of 
3.8  cm  and  subtend  an  angle  of  90°  as  shown  in  Fig.  2-14.  The  longi¬ 
tudinal  separation  of  adjacent  plate  edges  is  3.8  cm  with  a  center- 
to-center  spacing  of  7.6  cm.  The  azimuthal  angle  subtended  by  corres¬ 
ponding  coupling  plates  on  delay  line  1  and  delay  line  2  is  90°.  The 
edge  of  the  last  plate  is  located  approximately  12.7  cm  from  the  end- 
plate  of  the  experiment. 


Fig.  2-13  RF  Spectrum  of  Transmitter  Output 

First  peak  on  left  edge  is  the  zero  reference 
Log  vertical  scale 
2  MHz/division  horizontal  scale 
Signal  frequency  3.25  MHz 
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Chapter  3 

PLASMA  DIAGNOSTIC  SYSTEMS 
3.1  General  Description 

A  variety  of  diagnostic  systems  are  employed  to  investigate  plasma 
parameters  as  well  as  the  interaction  of  the  plasma  with  electrostatic 
waves  excited  by  the  slow  wave  structure.  The  diagnostic  svstems, 
which  will  be  discussed  in  this  chapter,  consist  of:  (1)  electrostatic 
(Langmuir)  probes,  (2)  dc  current  probes,  (3)  rf  probes,  and  (4)  elec¬ 
trostatic  energy  analyzer. 

Fig.  3-1  shows  the  basic  diagnostic  systems  used,  and  Fig.  3-2 
shows  the  locus  of  the  radial  and  azimuthal  probe  paths.  The  positions 
of  the  radially  and  azimuthally  movable  probes  are  accurately  determined 
by  mechanically  coupling  the  probe  assembly  to  a  10-turn,  1  kft  helipot. 
The  wiper  of  the  helipot  gives  a  voltage  proportional  to  the  angle  of 
the  probe  and  this  voltage  is  used  to  drive  the  x-axis  of  an  x-y  re¬ 
corder  as  shown  in  Fig.  3-3. 


3.2  Electrostatic  Probes 


Electrostatic  probes  are  frequently  used  for  the  localized  measure¬ 
ment  of  electron  temperature,  electron  number  density,  and  plasma  poten¬ 
tial  in  low  pressure  discharges.  The  theories  of  electrostatic  probes 
18  19  20  21 

are  well  documented  ’  ’  ’  and  only  the  relevant  equations  will 

be  given  here. 


The  construction  of  a  cylindrical  electrostatic  probe  is  shown  in 
Fig.  3-4  and  the  electronic  svstem  used  for  measurements  of  the  probe 
characteristic  and  the  radial  profile  of  ion  saturation  current  is  shown 
in  Fig.  3-5.  Fig.  3-6  shows  a  typical  probe  characteristic.  The  elec¬ 
tron  current  can  be  extracted  from  the  total  probe  current  when  Ad>  <  0, 
where 


(3.1) 


Here,  4>  and  represent  the  probe  potential  and  the  plasma  potential, 
respectively.  When  A<i<0,  that  is,  in  the  potential  range  where  elec¬ 
trons  are  being  repelled,  the  electron  current  is  given  by 


Fig.  3-1  Basic  Diagnostic  Setup 


Azimuthal  Probe  Path 


.ocus  of  Radial  and  Azimuthal  Probe  Paths 
(All  Dimensions  in  Centimeters) 
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Fig.  3-3  Circuitry  for  Probe  Positioning  Assembly 


Torr  Seal 


Fig.  3-4  Construction  of  Langmuir  Probe 
(Dimensions  in  Centimeters) 


f 
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Rg.  3-5  Electronic  Circuit  for  Langmuir  Probe  Measurements 


Rg.  3-6  Typical  Langmuir  Probe  Characteristic 
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where  ng  is  the  electron  number  density  and  k^T^/e  is  the  electron 

temperature  in  eV.  A'  represents  the  electron  current  collection 

P 

area.  Since  the  electron  current  is  constrained  to  flow  along  the 
background  magnetic  field  Bq,  the  electron  collection  area  is  approx¬ 
imated  by  the  projection  of  the  probe  area  on  a  plane  perpendicular 
to  Bq  rather  than  by  the  real  probe  area.  If  d  is  the  probe  diameter 
and  i  is  the  probe  length,  the  electron  collection  area  A'  is  given  by 


A'  -  2d£ 
P 


From  Eq.  (3.2)  it  is  found  that 


(In  i  ) 
e 


-1 


(3.3) 


(3.4) 


A  plot  of  versus  <£  yields  a  straight  line  whose  slope  determines 

the  electron  temperature.  The  electron  number  density  can  then  be  de¬ 
termined  by  Eq.  (3.5). 
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The  magnitude  of  the  ion  saturation  current,  |l  |  is  proportional  to 

1/2  is 

n^k^T^)  and  the  radial  profile  of  this  quantity  in  conjunction  with 

probe  characteristics  at  various  radii,  yield  valuable  Information  about 

the  radial  inhomogeneity  of  the  plasma.  These  characteristics  are  also 

measured  at  various  axial  locations  in  the  plasma  to  obtain  information 

about  the  axial  inhomogeneity  of  the  plasma. 


3. 3  DC  Current  Probes 


To  obtain  some  information  about  the  relative  axial  electron 
flux  over  the  cross  section  of  the  plasma,  both  radial  and  azimuthal 
current  probes  are  used  near  the  endplate  of  the  experiment.  The 
construction  of  these  current  probes  are  shown  in  Fig.  3-7  and  the 
measurement  circuit  is  shown  in  Fig.  3-8.  The  normal  to  the  surface 
of  the  current  probes  are  oriented  parallel  to  Bq  so  the  axial  com¬ 
ponent  of  electron  flux  is  measured.  The  diameter  of  the  circular 
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current  probes  are  made  small  so  as  Co  increase  Che  special  resolution 

of  the  probes.  The  probe  diameter  is  0.635  cm.  The  local  electron 

current  density  can  be  approximated  by  J  *  I  /A  where  I  is  the 

ez  cp  cp  cp 

measured  current  through  the  external  50Q  current  sensing  resistor 
and  A  is  the  cross  sectional  area  of  the  probe.  In  general,  Jgz  is 
a  function  of  r  and  9,  however  the  azimuthal  position  of  the  radially 
movable  current  probe  is  not  independent  of  the  probe  angle  a.  Be¬ 
cause  of  this  it  is  not  possible  to  extract  the  radial  dependence  of 
at  a  fixed  azimuthal  position.  It  is  found  however,  that  the  0 
dependence  of  Jgz  is  nearly  constant  in  the  azimuthal  range  of  9  near 
the  coupler  plates  (~<  9  < -^  ,  <  8  <~ ■  )  and  hence  the  a.  dependence 

of  J  can  be  separated  from  the  8  dependence.  The  total  current  picked 
ez 

up  by  the  endplate  of  the  device  is  determined  by  the  integration  of 
the  local  current  density  over  the  cross  section  of  the  plasma 


e, total 
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where  r^  is  the  plasma  radius. 

Prom  the  geometry  indicated  in  Fig.  3-2  the  radial  position  of  the  probe 
is  related  to  the  probe  angle  a  by. 


r  »  2R  sin  ■=■ 
o  2 


where  Rq  is  the  radius  of  the  arc  swept  out  by  the  probe.  Eq.  (3.7) 
is  also  used  to  relate  r  and  a  in  the  Langmuir  probe  and  rf  probe  meas¬ 
urements.  A  change  of  variables  from  (r,9)  to  (a, 6)  can  be  applied  to 
Eq.  (3.6)  yielding 


2  it  max 


e, total 


2Vsi”f)  Jez<a,9)  ifeU 


L  3(r,9) 

Where  TfaTe)  is  the  Jacobian  of  the  transformation  (r,9)  — »  (a, 9) .  The 
Jacobian  is  defined  by 
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Hence , 


3(r  ,9) 
3(o,8) 
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where  a  is  defined  by  Eq .  (3.11), 


~{k) 


(3.11) 


The  externally  measured  current  received  by  the  current  probe  is  in¬ 
fluenced  by  the  characteristics  of  the  plasma,  which  is  in  turn  influ¬ 
enced  by  the  structure  of  the  electromagnetic  fields  exciting  the  plasma. 
These  measurements ,  in  conjunction  with  the  electrostatic  energy  analyzer 
measurements  give  a  variety  of  information  regarding  the  axial  electron 
flux  and  the  axial  electron  energy  distribution  as  a  function  of  the  rf 
power  exciting  the  slow  wave  structure.  Experimental  results  and  their 
Interpretation  will  be  discussed  in  chapters  4  and  5  respectively. 

3.4  RF  Probes 

RF  prohes  have  been  extensively  used  in  plasma  physics  experiments 
to  investigate  amplitude  and  phase  characteristics  of  waves  in  plasma. 

9 ? 

For  example,  rf  probes  have  been  used  by  Briggs  and  Parker  “  in  the 

study  of  the  coupling  of  rf  energy  to  lower  hybrid  waves,  Fisher  and 
23 

Gould  in  the  investigation  of  the  formation  of  resonance  cones  in 
plasmas,  Bellan  and  Porkolab  in  the  study  of  the  excitation  of  lower 
hybrid  waves  using  a  slow  wave  structure,  and  Hooke  and  Bernabei  ^  in 
the  observation  of  propagating  lower-hybrid  waves. 

The  rf  probe  is  constructed  as  shown  in  Fig.  3-9.  As  can  be  seen 
in  the  figure  the  probe  is  doubly  shielded  (triaxial)  with  the  outer 
stainless  steel  tube  providing  the  vacuum  seal.  The  stainless  steel  tube 
is  inserted  through  a  probe  port  in  the  endplate  of  the  device  which 
allows  probe  motion  along  z  and  can  rotate  the  probe  in  a.  The  coaxial 


probe  tip  is  sealed  in  an  alumina  (A^O^)  insulating  shield  which  pre¬ 
vents  Che  probe  from  drawing  current  from  the  plasma.  The  amplitude 
and  phase  of  the  rf  potential  can  be  measured  using  the  electronic  sys¬ 
tem  illustrated  in  Tig.  3-10.  The  probe  output  is  connected  to  an  HP 
model  8405A  vector  voltmeter  and  a  50fl  termination  located  at  the  vec¬ 
tor  voltmeter  signal  input.  A  reference  signal  is  applied  through  a 
Tektronix  P6012,  10X  probe  from  one  of  the  slow  wave  structure  coupler 
plates  (see  Pig.  3-10)  to  the  reference  channel  on  the  vector  voltmeter. 
The  veccor  voltmeter  measures  the  amplitude  of  the  probe  signal  and  the 
phase  difference  between  the  probe  and  reference  at  the  frequency  of 
the  reference  signal.  Two  dc  output  voltages  are  available  from  the 
vector  voltmeter;  one  is  proportional  to  the  amplitude  of  the  probe 
signal  and  the  other  is  proportional  to  the  phase  difference  between 
the  probe  and  the  reference.  Each  of  these  voltages  is  applied  to  the 
y-axis  of  two  x-y  recorders  and  the  x-axes  are  driven  bv  the  probe  posi¬ 
tion  assembly  described  in  section  3.1.  Using  this  method  one  can 
accurately  map  the  amplitude  and  phase  of  the  rf  potential  as  a  function 
of  a  in  the  plasma.  These  measurements  are  then  made  at  various  axial 
locations  to  determine  amplitude  and  phase  characteristics  along  the 
plasma  column.  The  axial  position  is  measured  simply  bv  using  an  accu¬ 
rate  caliper. 

Care  must  be  taken  in  the  interpretation  of  the  phase  characteris¬ 
tics  since  the  vector  voltmeter  only  measures  phase  differences  between 
-180°  and  +180®.  If  the  phase  varies  by  more  than  360°  the  phase  out¬ 
put  of  the  vector  voltmeter  may  change  discontinuouslv  and  this  must  be 
distinguished  from  a  rapid  physical  change  in  phase.  Fortunately,  the 
phase  cha  ;es  are  generally  less  than  360°  across  the  plasma  column  and 
this  problem  is  not  experienced.  The  vector  voltmeter  also  possesses  a 
phase  offset  control  so  it  is  possible  to  bias  the  phase  output  at  a  con¬ 
venient  point  so  as  to  avoid  the  'ends'  of  the  phase  range  on  the  vector 
voltmeter  where  the  phase  output  may  change  discontinuouslv.  When  making 
phase  measurements  it  is  important  to  re<'r,.-v.  the  value  of  the  phase  off¬ 
set  to  make  comparisons  between  measurerae;* . 
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Fig.  3-9  Construction  of  RF  Probe 
(Dimensions  in  Millimeters) 
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Fig.  3-10  Electronic  Circuit  for  Amplitude  and  Phase  Measurements  of  RF  Potential 


32 


From  the  measured  phase  variation  as  a  function  of  a,  'jj(a) ,  it 
is  possible  to  determine  the  local  radial  wave  number  kf(r).  The 
radial  wave  number  can  be  readily  deduced  assuming  a  small  phase  var¬ 
iation  in  the  azimuthal  direction  near  the  coupler  plates  as  discussed 
in  conjunction  with  the  current  probes  in  section  3.3.  The  radial 
wave  number  is  determined  from  the  phase  information  by  Eq.  (3.12) 

k  (r)  =  ~  *(r)  (3.12) 

r  ar 


where  r  =  r(a)  is  defined  in  Eq.  (3.7).  k^_(r)  can  be  written  as, 


From  Eq.  (3.7)  we  find 


k  (r)  = 
r  da  dr 


dr  D  a 
—  =■  R  cos-r 
da  o  2 


(3.13) 


Hence, 


kr(r) 


_  a  \da 
R  cos-r 
o  2 


dd>  ,  , 

-7Z  .  a  *  ±  * 


(3.14) 


3.5  Retarding  Field  Electrostatic  Energy  Analyzer 

The  retarding  field  electrostatic  energy  analyzer  has  been  used  in 
plasma  physics  experiments  fo-  the  measurement  of  charged  particle  dis¬ 
tribution  functions  and  distribution  function  modification  due  to  the 

0  A 

presence  of  waves  in  the  piasma.  For  example,  Mau  "  has  studied  the 

theory  of  energy  analyzer  operation  and  the  modification  of  the  electron 

25  ’6 

distribution  due  to  damped  whistler  waves,  Andersen,  et.al.  *  has 
studied  the  ion  distribution  and  wave-particle  interaction  in  a  Q-machine, 
and  Guillemot,  et.al.  reported  the  deformation  of  the  electron  dis¬ 

tribution  due  to  damped  electrostatic  waves.  Design  considerations  for 
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energy  analyzers  were  reported  by  Simpson  . 

The  energy  analyzer  used  in  our  experiments  has  four  electrodes  as 
shown  in  Fig.  3-11.  The  entrance  grid,  shields  the  plasma  from  the 
electric  fields  within  the  analyzer  and  can  be  biased  at  the  plasma  poten¬ 
tial  p  assuming  this  does  not  disturb  the  operation  of  the  plasma.  Since 
this  grid  is  electrically  connected  to  the  endplate  of  the  experiment  and 

<h  is  a  function  of  radius  it  is  imoossible  to  bias  this  grid  at  the  local 
P 


plasma  potential  without  affecting  the  remaining  regions  of  the  plasma. 

In  our  case  it  is  connected  through  a  50P.  resistor  to  the  ECRH  coupler 
assembly  which  serves  as  the  reference  electrode.  Grid  2  is  the  dis¬ 
criminator  grid  and  imposes  an  axial  electric  field  to  retard  incoming 
charged  particles.  It  was  biased  at  negative  voltages  in  our  experi¬ 
ment  since  we  wanted  to  analyze  the  axial  energy  of  electrons.  If  $ 
is  the  local  plasma  potential  and  is  the  discriminator  potential 
only  electrons  with  energies,  <5>e(p^  -  Pj)  can  surpass  the  potential 
barrier  between  the  plasma  and  discriminator  grid  assuming  the  po- 
tential  on  the  entrance  grid  is  small.  It  is  experimentally  determined 
that  pg  is  sufficiently  small  ( — lOmV)  so  that  it  can  be  neglected. 

Grid  3  is  the  repeller  grid  and  is  biased  at  a  sufficiently  large  posi¬ 
tive  potential  so  as  to  repel  all  ions.  The  last  element  in  the  analyzer 
is  the  collector  electrode  which  collects  the  electrons  that  have  sur¬ 
passed  the  retarding  potential  barrier  of  the  discriminator  grid. 

The  construction  of  the  electrostatic  energy  analyzer  is  shown  in 
Fig.  3-12.  The  entrance  aperture  of  the  analyzer  is  located  at  a  radius 
of  3.81  cm  and  at  an  azimuthal  position  of  9=*ir/2.  i.e.  in  the  region 
of  maximum  signals  as  indicated  by  the  radial  and  azimuthal  current  probe 
measurements.  The  entrance  aperture  diameter  of  the  analyzer  is  1.11  cm. 
The  wire  .aesh  grids  were  constructed  from  tungsten  screening  and  the 
copper  collector  electrode  was  covered  with  a  layer  of  carbon  to  help 
eliminate  secondary  emission  from  the  collector.  Pertinent  dimensions 
of  the  energy  analyzer  are  indicated  in  Fig.  3-12. 

An  important  figure  of  merit  for  the  energy  analyzer  is  the  energy 

17  79 

resolution  which  has  a  lower  limit  determined  by  the  lens  effect  . 

The  entrance  aperture  of  the  energy  analyzer  can  be  considered  as  a  di¬ 
verging  electrostatic  lens  whose  focal  length  is  given  by 


S2  -  hYl 


(3.15) 


where  and  E^  are  the  electric  fields  in  the  plasma  and  in  the  analyzer 
respectively.  The  incoming  particles  have  an  energy  e&z-  The  electro¬ 
static  lens  is  diverging  since  E7  <  E^.  A  lower  bound  on  the  resolution 
can  be  estimated  by  assuming  E^  *  0  and  E^^O.  The  particle  energy  is 
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Entrance  Discriminator  Repeller  Collector 

Grid  Grid  Qrld  Electrode 


(See  Fig.  3-13) 


Fig.  3-11  Grid  Configuration  of  Retarding  Field  Electrostatic  Energy  Analyzer 


Fig.  3-12  Construction  of  Electrostatic  Energy  Analyzer 
(Dimensions  in  Centimeters) 
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»  »E  h  where  h  is  the  distance  between  the  entrance  grid  aperture  and 
the  discriminator  grid.  Hence  Eq.  (3.15)  reduces  to 

f  »-4h  (3.15A) 
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The  divergence  angle  ’  is  given  by 


,  2a 
sin  9  , 


(D/2) 

2  2 

(D/2)z  +  r 


(3.16) 


and  using  Eq.  (3.15A)  the  maximum  divergence  angle  is  given  by 


sin2  9 


„  (D/2)4 

d,max  (d/2)2  + 16h2 


(3.17) 


where  D  is  the  diameter  of  the  entrance  aperture.  It  is  now  assumed 
that  an  incoming  particle  possesses  both  axial  and  transverse  com¬ 
ponents  of  velocity,  v^  and  v^  respectively.  The  magnetic  field  does 
not  complicate  the  situation  since  both  electron  and  ion  gyro  radii 
are  much  less  chan  the  diameter  of  the  entrance  aperture.  The  trans¬ 
verse  energy  and  che  axial  energy  of  the  particle  is  related  by 

*1  2 

-=—  =•  tan  9 

^z 

where  9  Is  the  angle  between  the  velocity  vector  and  the  z-axis.  Since 
only  axial  momentum  is  sampled  by  the  analyzer,  the  energy  resolution 
is  given  by 


V1  .  2 
-=-  =  — ~  =  sin  9 

®  vZ 


(3.18) 


Hence  an  estimate  of  the  resolution  of  the  energy  analyzer  is  given 


(3.19) 


1+  (8h/D)  ‘ 


For  our  analyzer,  D*l.ll  cm  and  h*1.43  cm  which  gives  a  resolution 
A5/5  =  0.009. 

The  energy  analyzer  is  designed  to  sample  electrons  with  a  posi¬ 
tive  axial  velocity  component  (directed  away  from  the  ECRH  plasma  source 
The  electrons  possess  an  axial  distribution  function  fg(v  )  an<^  i-s 
desired  to  find  a  relationship  between  the  collector  current  1^,  the  dis 
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criminator  potential  and  the  distribution  function  f  (v  ) •  The 

a  q  z 

collector  current  is  determined  by  the  number  of  electrons  that  sur¬ 
passed  the  potential  barrier  established  by  the  discriminator  grid. 
Hence , 

Ic(*d)  -  “ATfe  |  vzfe(vz)dvz  (3.20) 

v 

c 

where  vz is  the  axial  electron  velocity  associated  with  an  energy 


-e$  »(l/2>n vz  ,  A  is  the  cross  sectional  area  of  the  entrance  aperture, 

T  is  the  overall  geometric  transmission  coefficient  of  the  grids,  and 
v^  is  the  cutoff  velocity  determined  by  the  discriminator  and  plasma  po¬ 
tentials. 


VO 

Z  o 


dv 


/-2e<A 

(  0 
\  m  / 
\  e  / 

(^) 


dd 


(3.21) 


(3.22) 


Substituting  Eqs.  (3.21)  and  (3.22)  into  Eq.  (3.20)  yields 

U*  'l/2 

o 


V*d>  *  -ATe 


\l/2 


,  2m  4> 

1  e  oi 


P  a 


I  <40 


-ATe 


c  a  m 


f  (v  (4>  ))d4> 
e  z  o  o 


(3.23) 


P  a 


Differentiating  yields, 


(v  (<t>  - 
z  p 


40) 


(3.24) 


Eq.  (3.24)  demonstrates  that  the  distribution  function  can  be  obtained 
by  differentiating  the  collector  current  with  respect  to  the  discriminator 
potential. 
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The  differentiation  of  Eq.  (3.24)  is  performed  electronically  as 
17,27 

shown  in  Fig.  3-13  ’  .  It  was  found  experimentally  that  a  repeller 

grid  potential,  ^R«100V  was  sufficient  to  repel  all  ions.  The  output  of 
the  collector  electrode  was  connected  through  a  shielded  coaxial  cable 
to  the  input  of  a  Princeton  Applied  Research  model  JB-5  lock-in  amplifier 
(input  impedance  of  50kft) .  The  discriminator  voltage  was  swept  by  the 
sawtooth  output  (0-  +  150V)  from  a  Tektronix  type  33  oscilloscope.  A 
-150V  power  supply  was  kept  in  series  with  this  sawtooth  signal  to  main¬ 
tain  the  discriminator  at  a  negative  potential  since  we  are  sampling  the 
electron  energies.  Also  a  small  amplitude,  high  frequency  reference 
potential,  was  applied  to  the  discriminator  grid  through  an  isola¬ 

tion  transformer.  As  can  be  seen  in  Fig.  3-13  the  reference  was  derived 
from  the  internal  reference  of  the  lock-in  amplifier. 

Reference  Signal:  $  ,(t)-v  cosw  ,t 

ret  o  i  _-r 

DC  +  Sawtooth  Signal:  4  (t)  =*  -V  +V  cosw  .  ,t 

s  o  o  s ,fund 


us,fundis  cJ,e  fundamental  frequency  component  of  the  sawtooth  signal. 
The  total  discriminator  potential  is  given  by 


Conditions  were  set  on  the  applied  signals  such  that 


V»v,u»  ,»u)  , 

o  o  ref  s,fund 

The  time  base  on  the  oscilloscope  was  typically  set  at  a  sweep  of  5 
sec/cm  and  the  reference  frequency  from  the  lock-in  amplifier  at  1  kHz. 
On  this  time  scale  the  voltage  sweep  from  the  sawtooth  generator  is  ex¬ 
tremely  slow  compared  to  the  time  scale  of  the  reference  signal.  The 
collector  current  is  written  as, 

I  (♦,)  *  I(*  +f  f(t))  (3.25) 

c  d  c  s  ref 

Since  V  >>  v  ,  we  expand  I  (4.)  about  4  retaining  terms  only  to  first 

0  0  C  Q  S 

order 


Fig.  3-13  Electronic  System  for  Measuring  Electron  Distribution  Function 
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The  collector  signal  is  mixed  with  the  reference  signal  in  the 
lock-in  amplifier  and  a  dc  component  is  produced  due  to  the  term 


‘ref"5 


s 


This  dc  output  is  filtered  through  a  low  pass  filter  with  a  time  con¬ 
stant  of  1  sec  (fast  compared  with  sweep  speed)  and  is  applied  to  the 
y-axis  of  an  x-y  recorder.  The  x-axis  is  driven  by  the  sawtooth  sweep 
signal  $  .  Using  this  technique  the  energy  distribution  function  of 
electrons  can  be  displayed  directly  on  the  x-y  recorder.  It  should  be 
noted  that  the  reference  signal  amplitude  v  must  satisfy  the  condition 

that  v  <<  V  ;  it  should  not  be  so  small  as  to  sample  the  small  scale 
o  o 

fluctuations  in  the  distribution  function.  The  entire  electronic  sys¬ 
tem  is  isolated  from  ground  through  an  isolation  transformer  which 
supplies  all  ac  power  to  the  electronic  equipment.  The  ECRH  coupler 
serves  as  the  common  reference  electrode. 


The  electron  velocity  distribution  is  normalized  as 
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<n> 

total 


j  fe(Vz)dVz 
—  00 


(3.27) 


Some  information  regarding  the  dynamics  of  the  electrons  can  be  de¬ 
termined  by  analyzing  the  moments  of  the  electron  distribution.  Care 
must  be  taken  in  the  interpretation  of  the  energy  analyzer  data  since 
the  moments  as  determined  from  the  experimental  distribution  are  not  the 
same  as  the  moments  obtained  from  the  real  distribution  since  the  energy 
analyzer  does  not  sample  electrons  with  negative  velocities.  A  typical 
experimental  electron  distribution  is  shown  in  Tig.  3-1^  where  $  is  the 

absolute  value  of  the  applied  dc  discriminator  potential  6  and  is  the 

s  p 

negative  of  the  plasma  potential.  The  moments  of  the  experimental  dis¬ 
tribution  g(<&)  are  given  in  Eqs.  (3.28)  to  (3.33)  i?,29 
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If  <n>  is  uncorrelated  with  <nv^>  and  <—  raenvz>  the  following 


can  also  be  written 
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and  the  mean  random  energy  <<5Z>  can  be  written 


(3.28) 


(3.29) 


(3.30) 

moments 


(3.31) 


(3.32) 
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Eqs.  (3.28)  to  (3-33)  are  evaluated  numerically  for  the  experimental  data. 

It  should  be  noted  that  if  there  is  a  drift  velocity  associated  with 
the  electron  distribution  and  the  plasma  potential  is  not  known  a  priori 
it  may  be  difficult  to  interpret  the  experimental  data  since  a  change  in 
plasma  potential  also  will  shift  the  experimental  distribution  function  as 
well  as  a  shift  in  energy  of  the  real  distribution  function.  However,  the 
combination  of  energy  analyzer  data,  current  probe  measurements,  and  electro¬ 
static  probe  measurements  will  help  in  the  understanding  of  the  electron 
dynamics  as  will  be  discussed  in  chapter  4. 
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Chapter  4 

EXPERIMENTAL  RESULTS 

4.1  Introduction 

To  understand  the  wave-plasma  interaction,  it  is  important  to  be 
familiar  with  the  operating  characteristics  of  the  plasma  device,  the 
unperturbed  nature  of  the  plasma,  and  the  free  space  characteristics 
of  the  slow  wave  structure.  By  doing  this,  it  is  possible  to  obtain 
a  better  understanding  of  how  the  rf  electrostatic  waves  interact  with 
the  plasma  and  modify  it.  In  this  chapter  a  discussion  of  the  back¬ 
ground  plasma  and  Che  rf  experiments  will  be  presented, and  in  Chapter  5 
an  interpretation  of  these  results  along  with  a  simple  theoretical  model 
will  be  presented. 

All  plasma  experiments  were  performed  using  Argon  at  an  operating 
-4 

pressure  of  2  x  10  Torr  and  the  axial  magnetic  field  profile  illus¬ 
trated  in  Fig.  2-4. 

4.2  Background  Plasma 

The  incident  ECRH  microwave  power  determines  the  initial  state  of 
the  plasma  before  low  frequency  electrostatic  waves  are  excited  by  the 
slow  wave  structure.  Electrons  absorb  microwave  energy  through  ECRH 
and  the  hot  electron  fluid  expands  out  of  the  coupler  region  into  the 
drift  region  with  a  relatively  uniform  axial  magnetic  field.  The  cold 
ions  then  acquire  an  axial  drift  out  of  the  coupler  region  to  maintain 
quasineutrality.  As  the  microwave  power  is  increased,  the  ionization, 
and  hence  the  local  electron  density  increases  until  the  density  reaches 
a  level  such  that  the  local  electron  plasma  frequency  is  equal  to  the 
incident  microwave  frequency.  At  this  point  additional  microwave  energy 
is  shielded  from  the  plasma  and  is  reflected  from  the  microwave  coupler 
and  is  dissipated  in  the  isolator.  Hence,  the  maximum  achievable  electron 
density  is  determined  from  this  criterion.  The  electron  plasma  frequency 
is  determined  from  the  electron  density  by  Eq.  (4.1) 


For  the  microwave  frequency  of  2.45  GHz  used  in  the  experiment  the  upper 
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limit  on  electron  density  is  7.43  x  10  cm  .  Fig.  4-1  shows  a  graph 
of  the  reflected  and  transmitted  microwave  power  for  the  Lisitano  coil 
versus  incident  microwave  power.  It  can  be  seen  in  Fig.  4-1  that  ap¬ 
proximately  25  watts  is  required  to  initiate  a  discharge.  It  should 
also  be  noted  that  the  reflection  coefficient  is  unity  if  the  plasma 
is  not  present  due  to  the  short  circuit  termination  on  the  Lisitano 
coil.  For  most  experiments  the  system  is  operated  at  an  incident  micro- 
wave  power  of  490  watts  and  a  transmitted  power  of  121  watts. 

Fig. 4-2  defines  the  axial  reference  position,  z  **  0,in  reference 
to  the  slow  wave  structure.  It  should  be  noted  that  the  z-scale  used 
here  is  not  the  same  as  the  z-scale  used  in  some  of  the  figures  of 
chapter  2.  This  axial  reference  is  important  to  define  since  the  plasma 
is  inhomogeneous  in  both  the  axial  and  the  radial  directions  and  both 
inhomogeneities  affect  the  coupling  of  the  rf  energy  from  the  slow  wave 
structure  to  the  plasma.  Fig.  4-3  shows  the  radial  profile  of  the  elec¬ 
tron  density  at  an  axial  position  of  z  *  -3.81  cm.  Figs.  4-4  and  4-5 

are  semi-log  plots  of  the  electron  density  and  electron  temperature 
2 

versus  r  at  z  =  -3.81  cm.  A  radial  profile  of  the  relative  magnitude 
of  the  ion  saturation  current  is  shown  in  Fig. 4-6.  Figs.  4-3  and  4-6 
indicate  that  the  radial  profile  of  ion  saturation  current  is  approxi¬ 
mately  the  same  as  the  radial  profile  of  electron  density  indicating 
that  the  electron  temperature  varies  more  slowly  across  the  plasma 
than  does  the  electron  density.  This  result  is  supported  by  Figs.  4-4 

and  4-5  which  indicates  the  radial  scale  length  of  electron  temperature 
variation  is  longer  than  the  radial  scale  length  of  the  electron  density 
variation.  Fig.  4-7  shows  the  magnitude  of  the  ion  saturation  current 
versus  the  radial  probe  angle  a  with  the  axial  position  z  as  a  parameter. 
From  Fig.  4-7  and  various  probe  characteristics  at  different  axial  posi¬ 
tions  in  the  plasma, the  axial  profile  of  electron  density  and  electron 
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Fig.  4-1  Reflected  and  Transmitted  Microwave  Power  vs.  incident  Power  for  Lisitano  Coil 


(Dimensions  in  Centimeters) 


Fig.  4-3  Radial  Profile  of  Electron  Density  at  z  =  -3.81  cm 


(cm“3) 


Fig.  4-7  Magnitude  of  Ion  Saturation  Current  vs.  a  with  z  as  a 
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temperature  can  be  determined.  Figs.  4-8  and  4-9  show  the  axial  pro¬ 
files  of  electron  density  and  electron  temperature  at  r *  0  respectively. 

The  Langmuir  probe  measurements  also  indicate  that  the  background  plasma 
is  azimuthally  symmetric.  Fig.  4-10  shows  the  magnitude  of  the  ion  sat¬ 
uration  current  versus  the  radial  probe  angle  a  with  the  transmitted  ECRH 
power  as  a  parameter.  The  effect  of  the  reduction  of  ECRH  power  is  to 
reduce  the  ion  saturation  current  and  hence  the  electron  density  profile 
as  well  as  the  electron  temperature  profile.  It  is  found  experimentally, 
however,  that  the  electron  temperature  decreases  at  a  slower  rate  than 
does  the  electron  density  with  ECRH  power  and  hence  the  curves  in  Fig.  4-10 
are  indicative  of  the  electron  density  profile  in  the  plasma.  Based  on 
these  Langmuir  probe  results,  it  is  possible  to  construct  a  physical 
model  of  the  background  plasma  electron  density. 

The  spatial  variation  of  electron  density  is  modeled  by 


n  (r,z)  =  n  e_r  (1-f)  (4.2) 

where  R  and  L  are  the  radial  and  axial  scale  lengths  determined  from 
Figs.  4-4  and  4-8  respectively,  n^  is  the  electron  density  at  the 
position  (r,z)  «  (0,0).  These  three  parameters  are  given  by: 

n  -  4.82  x  1015m~3 
o 

R  "  0.03ra 
L  *  0 . 39m 

Radial  and  azimuthal  current  probes  were  used  to  measure  the 
localized  current  density  in  the  plasma.  The  electron  current  density 
Jez  versus  the  radial  probe  angle  a  is  shown  in  Fig.  4-11  and  the  ex¬ 
perimental  results  also  indicated  that  the  localized  current  density 
was  azimuthally  symmetric  as  determined  from  the  azimuthal  current 
probe.  The  total  current  drawn  from  the  endplate  of  the  device  is 

approximately  -2  mA  and  a  numerical  integration  of  J  versus  a  yields 

ez 

approximately  -1.8  mA  comparing  favorably  with  the  measured  value. 


f 


—  


Fig.  4-8  Axial  Profile  of  Electr 


10  Magnitude  of  Ion  Saturation  Current  vs.  a  with  ECRH  Power  as  a  Parameter 
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The  electrostatic  energy  analyzer  is  located  at  a  radial  position 
of  3.81  cm  and  an  azimuthal  position  of  9  =  as  discussed  in  section 
3.5.  The  endplate  of  the  experiment  and  hence  the  entrance  aperture 
of  the  energy  analyzer  is  located  at  an  axial  position  of  z  *  29.8cm 
(z  *  0,  defined  in  Fig.  4-2).  The  axial  electron  energy  distribution 
is  shown  in  Figs.  4-12  and  4-13.  From  Fig.  4-13  the  axial  electron 
temperature  is  determined  to  be  4.31eV  which  compares  favorably  with 
the  Langmuir  probe  results.  The  maximum  value  of  the  electron  energy 
distribution  can  be  normalized  by  using  the  measured  electron  density 
as  determined  by  the  Langmuir  probe  and  a  numerical  integration  of 
Eq.  (3.28). 


4.3  Typical  Plasma  Parameters 


Table  4.1  is  a  summary  of  typical  plasma  parameters  at  (r,z)  = 
(0,0).  Coulomb  collision  frequencies  are  determined  from  Spitzer's 
formulas 


e  n  in  A 
e 


ee 


16Tr1/2eV/2(k  T  )3/2 
o  e  Be 


(4.3) 


where 
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re 


ii 


ee 
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12r(e  k  T  /e2)3/2 
o  B  e 


1/2 


(4.4) 

(4.5) 

(4.6) 

(4.7) 
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For  a  neutral  pressure  of  2  x  10  Torr  the  neutral  density  is 
The  collision  cross  section  for  an  electron-neutral 
Argon  collision  for  an  electron  temperature  at  (r,z)  ■  (0,0)  of  5.5eV 


1  8 

6.2  x  10  m 


is  9  x  10  m.  The  electron  thermal  speed  at  this  temperature  is 
approximately  106m/s.  The  electron-neutral  collision  frequency  is 


given  by 


v  »no  v 
en  n  ens  te 


Similarly  the  ion-neutral  collision  frequency  is  given  by 


v.  =  n  o.  v  . 
in  n  ms  ti 


(4.9) 


where  o . 


-19  2 
4  x  10  m. 


Additional  interaction  cross  sections  are  given  below 

-  Ionization  of  Argon  atoms  by  an  electron: 

1n-21  2 
a  .  -  10  m 
em 

9 

-  Excitation  of  “P  levels  of  Argon  atoms  by  an  electron: 

s 

ln-21  2 
a  '  10  m 
enx 

-  Charge  exchange  collisions  between  Argon  atoms  and  an  ion: 

-19  2 

a.  ,  '  4  x  10  m 
inch 


0 
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Table  4-1  Typical  Background  Plasma  Parameters  at  (r,z)  °  (0,0) 


Quantity  Svmbol 

Mass  of  Argon  Ion  nu 

Magnetic  Field  B 

o 

Neutral  Pressure  P 

n 

Neutral  Density  n 

n 

Electron,  Ion  Density  n 

6,1 

Radian  Electron  Plasma  Frequency  jj 

pe 

Electron  Plasma  Frequencv  f 

-  pe 

Radian  Ion  Plasma  Frequency 

Ion  Plasma  Frequency  f  ^ 

Radian  Electron  Cyclotron  Frequency 
Electron  Cvclocron  Frequencv  f 

e 

Radian  Ion  Cyclotron  Frequency  12^ 

Ion  Cyclotron  Frequency 

Electron  Temperature  'tBTe 

Ion  Temperature  kBTi 

Electron  Thermal  Speed  v 

Ion  Thermal  Speed  v  ^ 

Radial  Density  Profile  Scale  Length  R 

Axial  Density  Profile  Scale  Length  L 

Debye  Length 

Debve  Wave  Number  k, 

d 

Electron  Gyro  Radius  p 

e 

Ion  Gyro  Radius  p. 

Electron-Electron  Collision  Frequency 
Electron-Ion  Collision  Frequency  v  ^ 

Ion-Electron  Collision  Frequency  v ^ 

Ion-Ion  Collision  Frequency 
Electron-Neutral  Scattering  Rate  Vens 

Electron-Neutral  Ionization  Rate  Veni 

Electron-Neutral  Excitation  Rate  Ven> 

Ion-Neutral  Scattering  Rate  V in= 

Ion-Neutral  Charge  Exchange  Rate  v inc 

Radian  Signal  Frequency  w 

Signal  Frequency  f 


Value 


m. 

i 

6.63  x  10“26kg 

B 

850  Gauss 

0 

-4 

Pn 

n 

2  x  10  Torr 

,  .  . a18  -3 

6.2  x  10  m 

n 

, .15  -3 

n 

4.62  x  10  m 

e,i 

jj 

3.92  x  109  rad/i 

pe 

f 

624  MHz 

pe 

7 

a  . 

1.45  x  10  rad/< 

Pi 

f  . 

2.3  MHz 

Pi 

a 

1.50  x  1010  rad, 

e 

f 

2.39  GHz 

e 

ai 

2.05  x  109  rad/; 

fi 

32.6  kHz 

¥e/e 

5.5  eV 

k_T./e 

-  .5  eV 

a  i 

V 

10°  m/s 

tfc  -j 

V  . 

1.1  x  10  m/s 

tl 

R 

.03  m 

L 

.39  m 

XJ 

1.78  x  10"4  m 

d 

„3  -1 

k , 

35.2  x  10  m 

d 

-5 

p 

6.67  x  10  m 

e 

-3 

o  . 

5.37  x  10  m 

1 

V 

27  kHz 

ee 

Jei 

9.53  kHz 

v  . 

.13  Hz 

le 

uu 

3.64  kHz 

V 

558  kHz 

ens 

V 

6.2  kHz 

eni 

V 

6.2  kHz 

enx 

V 

2.8  kHz 

ms 

\i 

2.8  kHz 

inch  7 

a) 

2.04  x  10  rad / 

3.25  MHz 
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4.4  Free  Space  Characteristics  of  the  Slow  Wave  Structure 


The  physical  construction  of  the  slow  wave  structure  was  described 

in  section  2.6.  In  this  section  the  rf  potential  associated  with  the 

slow  wave  structure  in  free  space  will  be  described.  The  fundamental 

wave  number  associated  with  the  slow  wave  structure,  k  is  determined 

20 

by  Eq.4.10 


k  =  ^ 
20  Az 


(4.10) 


where  M >  is  the  phase  delay  per  section,  and  Az  is  the  distance  between 

adjacent  sections  of  the  slow  wave  structure.  For  the  experiments  de- 

scribed  here  —  and  Az  =  7.6  cm  yielding  !^ZQi=  *21  cm  .  Fig. 4-14 

shows  the  axial  phase  delay  as  measured  by  the  rf  probe  (see  Section  3.4) 
ik  z  -ik  z 

for  both  e  20  and  e  20  excitations  at  a  radial  position  located 
a  small  distance  away  from  the  coupler  plates.  At  a  radial  position 
corresponding  to  the  coupler  plates  it  would  be  found  that  over  the 
length  of  the  coupler  plate  there  would  be  no  phase  change  and  hence 
the  phase  change  along  the  axis  would  not  be  as  smooth.  This  is  due 
to  the  high  k^  content  in  the  spatial  Fourier  transform  of  the  slow 
wave  structure. 


Fig.  4-15  shows  the  relative  magnitude  of  the  rf  potential  versus 
the  radial  probe  angle  a  at  various  axial  locations  for  an  rf  power  of 
100  watts  in  the  delay  line  corresponding  to  the  coupler  plates  at  9  * 
Since  the  slow  wave  structure  is  a  balanced  system,  the  magnitude  of  the 
rf  potential  is  an  even  function  of  a.  The  radial  wave  number  is  imagin¬ 
ary  since  the  axial  wavelength  is  much  less  -''an  the  free  space  wave¬ 
length;  i.e.  —  <<  k  and  the  signal  is  evanescent  in  the  radial  dir- 

C  20 

ection.  Fig.  4-16  shows  the  a-profile  of  phase  and  it  is  seen  that  the 
phase  is  constant  except  for  a  nearly  discontinuous  180°  change  in  phase 
at  a  *  0.  The  constant  phase  fora^Ois  indicative  of  the  evanescence 
and  the  180°  phase  change  is  due  to  the  fact  that  the  structure  is 
balanced. 


N>|  =1 


Peah  Probe  Signal  (mV) 
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Prf  =  100  watts 


z  =  5.72cm 


z  =  5.40cm 


z  =3.8 1cm 
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Fig.  4-15  Relative  Magnitude  of  RF  Potential  vs.  a 
with  z  as  a  Parameter  in  Free  Space 
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4.5  Experimental  Results  on  the  Interaction  of  the  Plasma  with 
Unidirectional  R F  Waves  Excited  by  the  Slow  Wave  Structure 

This  section  is  devoted  to  a  detailed  compilation  of  experimental 
results  involving  the  excitation  of  the  plasma  by  rf  fields  from  the 
slow  wave  structure.  Unidirectional  slow  waves  propagating  on  the  ex¬ 
citing  structure  can  be  set  up  with  k  >  0  or  k  <  0  by  interchanging 
the  power  feed  points  and  the  terminating  50ft  loads  on  the  delay  lines  as 
discussed  in  section  2.6.  Experimental  results  are  given  for  both  k2Q>  0  and 
k  <  0  excitations.  Results  on  the  following  topics  are  included  in 
this  section;  (1)  the  radial  and  axial  profiles  of  electron  density, 
electron  temperature,  and  ion  saturation  current  for  different  rf  powers, 

(2)  the  radial  profile,  azimuthal  profile  and  rf  power  dependence  of 
axial  electron  flux,  (3)  the  affect  of  rf  power  on  the  electron  energy 
distribution  function,  and  (4)  the  radial  and  axial  profiles  of  the 
amplitude  and  phase  of  the  rf  potential  as  a  function  of  rf  power. 

The  results  of  the  rf  modification  of  the  electron  density  and 

ion  saturation  current  profiles  are  shown  in  Figs.  4-17  to  4-19. 

2 

Fig. 4-17  shows  the  electron  density  versus  r  and  Figs.  4-18  and  4-19 

show  the  ion  saturation  current  versus  o  for  different  rf  powers,  » 

as  well  as  for  k  >0  and  k  <  0. 

zo  zo 

Even  though  the  electron  density  profiles  are  not  Gaussian  for 

P  ^  f  0,  it  is  seen  from  Fig.  4-17  that  it  is  possible  to  assume  a 

local  Gaussian  profile  for  r  >_  2.5  cm  with  a  local  radial  scale  length 

equal  to  the  radial  scale  length  of  the  unperturbed  plasma.  It  should 

be  noted  that  the  absolute  electron  density  decreases  as  the  rf  power 

is  increased.  It  is  also  found  that  for  a  given  rf  power  the  density 

profile  for  a  k  <0  excitation  is  smaller  than  for  a  k  >0  excita- 
zo  zo 

tion.  These  results  are  corroborated  by  the  a-profiles  of  ion  satura¬ 
tion  current.  Experimental  results  also  indicate  that  the  axial  elec¬ 
tron  density  gradient  is  not  affected  along  the  slow  wave  structure. 

To  understand  the  wave-particle  dynamics,  it  is  important  to  be  familiar 
with  these  modified  density  profiles. 


Prf  =  6.25  watts 


Magnitude  of  Ion  Saturation  Current  vs.  a  for  Prf  =  6.25  Watts 


Magnitude  of  Ion  Saturation  Current  vs.  a  for  Prf  =  25  Watts 


Radial  and  azimuthal  dc  current  probes  are  used  to  measure  the 

profiles  of  axial  electron  flux  in  the  plasma  as  described  in  section 

3.3.  These  measurements  are  made  as  a  function  of  rf  power  and  the 

direction  of  excitation  of  the  slow  wave  structure,  a-profiles  of 

electron  flux  for  various  rf  powers  are  shown  in  Figs.  4-20  and  4-21 

for  k  ■  .21  cm  ^  and  k  ■  -.21  cm  ^  respectively.  The  azimuthal 
zo  zo 

electron  flux  at  a.  radius  of  3.8  cm  is  shown  in  rig.  4-22.  Fig.  4-22 

applies  for  both  k  >0  and  k  <0  excitations. 
r  zo  zo 

As  can  be  seen  in  Figs.  4-20  and  4-21  there  are  two  radially 

localized  regions  of  axial  electron  flux.  Plots  of  the  peak  electron 

current  density  for  peaks  A  and  B  (see  Figs.  4-20  and  4-21)  as  a 

function  of  rf  power  for  k  >0  and  k  <0  are  shown  in  Figs.  4-23 

and  4-24.  Fig.  4-25  shows  the  radial  location  of  peak  A  of  the  electron 

flux  as  a  function  of  rf  power  for  +  k  excitations. 

—  zo 

If  the  microwave  power  used  for  ECRH  is  reduced,  the  radial  pro¬ 
files  of  electron  flux  are  radically  altered.  Figs.  4-26  and  4-27  show 
the  effect  of  the  reduction  in  ECRH  power  on  the  a-profile  of  axial 

electron  current  density  for  k  >0  and  k  <0  respectively.  From 

zo  zo 

Fig.  4-10  it  is  seen  that  the  ion  saturation  current  and  hence  the 
electron  density  profiles  are  reduced  as  ECRH  power  is  reduced  and 
from  Figs.  4-26  and  4-27,  it  is  seen  that  peak  A  moves  radially  inward, 
and  peak  B  remains  in  approximately  the  same  radial  position.  Peak  A 
moves  inward  in  such  a  fashion  that  the  location  of  this  peak  always 
corresponds  to  the  same  value  of  electron  density.  The  electron  den- 

3  n 

sity  at  this  point  is  n^  -  6  x  10  cm  .  Peak  B  remains  approximately 
in  the  same  position  regardless  of  the  power  level  used  to  excite  the 
slow  wave  structure.  The  location  of  peak  B  also  is  not  affected  by 
the  electron  density  profile.  These  results  indicate  peak  A  is  due 
to  the  excitation  of  a  plasma  mode  by  the  slow  wave  structure  while 
peak  B  is  probably  due  to  the  large  axial  electric  fields  located  near 
the  coupler  plates.  In  this  report  we  will  be  mainly  concerned 
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Profiles  of  Axial  Electron  Current  Density  for  kzo=.21cm 


Profiles  of  Axial  Electron  Current  Density  for  kzo=-21cm 


Fig.  4-22  Azimuthal  Profile  of  Electron  Current  Density  at  r=  3.81cm 


Max.  Value  of  Electron  Current  Density  for  Peak  A  (A/m2  ) 
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Fig.  4-23  Max.  Value  of  Electron  Current  Density  for  Peak  A  vs.  RF  Power 


.  Value  of  Electron  Current  Density  for  Peak  B  (A/m2) 


r  (cm) 
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Fig.  4-25  Radial  Location  of  Peak  A  of  Electron 
Flux  as  a  Function  of  RF  Power 


Electron  Current  Density  vs.  a 
vith  ECRH  Power  as  a  Parameter 
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Fig.  4-27  Axial  Electron  Current  Density  vs.  a 
for  kZo  =  -.21  cm-1  with  ECRH  Power  as  a  Parameter 
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with  the  electron  flux  associated  with  peak  A.  From  Fig.  4-23  it  is 

seen  that  for  rf  power  levels  less  than  approximately  40  watts,  the 

electron  current  density  increases  linearly  with  applied  rf  power 

2  2 

with  slopes  of  0.028  A/m  -watt  and  0.014  A/m  -watt  for  k  >0  and 

zo 

k  <  0  respectively.  For  power  levels  greater  than  approximately 
100  watts  the  electron  current  density  saturates  for  a  kzo  <  0  exci¬ 
tation.  These  measurements  are  corroborated  with  measurements  using 
the  electrostatic  energy  analyzer. 

Additional  information  on  rf  driven  axial  electron  currents  is 

obtained  using  the  electrostatic  energy  analyzer.  A  superposition  of 

typical  electron  energy  distributions  with  rf  power  as  a  parameter 

are  shown  in  Figs.  4-28  and  4-29  for  kzQ  >  0  and  <  0  respectively. 

It  is  evident  that  for  both  k  >0  and  k  <0  that  there  is  a  shift 

zo  zo 

in  the  location  of  the  peak  of  the  distribution  function  to  higher 
energies  as  the  rf  power  is  increased.  Fig.  4-30  shows  the  shift  in 
the  location  of  the  peak  of  the  energy  distribution  function  versus 
rf  power  and  Fig. 4-31  shows  the  change  in  velocity  corresponding  to 
the  energy  peak  versus  rf  power.  If  in  g($)  is  plotted  versus  $  for 
4>  greater  than  the  energy  corresponding  to  the  peak  of  the  distribution 
a  straight  line  results  whose  slope  determines  the  effective  axial  elec¬ 
tron  temperature,  T  .  Fig. 4-32  shows  Tez  versus  the  applied  peak  rf 
voltage  on  the  coupler  plates.  The  relative  value  of  axial  electron 
flux  <nvz>  and  the  average  electron  velocity  <yz>  as  determined  from 
Eqs.  3.29  and  3.31  are  plotted  versus  rf  power  in  Figs.  4-33  and  4-34 
respectively.  From  Fig. 4-33  it  is  seen  that  the  axial  electron  flux 
increases  linearly  with  rf  powers  less  than  40  watts  at  a  rate  twice 
as  large  for  k^^  >  0  than  for  kzQ  <  0  corroborating  the  dc  current 
probe  results  shown  in  Fig. 4-23.  It  is  also  observed  that  <nv?> 
saturates  for  Prf  *  100  watts  for  k  <0.  Fig. 4-35  is  a  plot  of  the 

zo  i  2 

average  kinetic  energy  of  an  electron  <  -j  m^v  >  as  determined  from 
Eq.3.28  versus  rf  power  for  +  excitations.  From  Fig. 4-29  one  ob¬ 
serves  Che  formation  of  a  plateau  in  the  distribution  function 


Fig.  4-28  Electron  Energy  Distribution  for  k^^lcm1  with  Prf  as  a  Parameter 
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Fig.  4-29  Electron  Energy  Distribution  for  kz0  =  -.21cm  with  Prf  as  a  Parameter 
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Fig.  4-31  Velocity  Corresponding  to  Peak  of 
Energy  Distribution  vs.  Prf 
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Fig.  4-33  Relative  Value  of  Axial  Electron  Flux  vs.  RF  Power 


(s/m  oO  <  A > 


for  power  levels  greacer  than  1/0  watts.  Fig. 4-36  is  a  typical 
semi-log  plot  of  g ( <*>)  versus  $  showing  the  development  of  the 
plateau  for  energies  greater  than  approximately  50eV. 

Up  to  this  point  we  have  concentrated  on  an  experimental 
description  of  the  electron  dynamics.  We  will  now  concentrate 
on  the  experimental  results  from  the  rf  probe  which  will  yield 
some  information  on  the  wave  dynamics  in  the  plasma.  The  rf 
measurements  consist  of  amplitude  and  phase  profiles  at  various 
axial  locations  in  the  plasma.  The  ot-profile  of  the  amplitude  of 
the  rf  potential  is  shown  in  Fig.  4-37.  The  free  space  amplitude 
profile  is  also  shown  for  comparison.  Fig. 4-37  indicates  a  shield¬ 
ing  of  the  rf  field  by  the  plasma.  This  result  was  observed  in  a 

13 

numerical  simulation  by  Decyk,  et  al  .  A  sharp  dip  in  the  potential 
magnitude  is  seen  at  an  angle  of  a  *  27.5°.  The  maximum  value  of  the 
peak  rf  amplitude  is  shown  in  Fig. 4-38  as  a  function  of  the  peak 
exciting  voltage  between  adjacent  coupler  plates.  Experimentally  it 
is  found  that  the  magnitude  of  the  rf  potential  decreases  in  the  axial 
direction  near  peak  A.  Figs.  4-39  and  4-40  show  typical  a-profiles 
of  phase  for  k  >0  and  k  <0  excitations  respectivelv  at  different 
axial  posit-ons  for  an  rf  power  of  6.25  watts.  Similar  profiles  are 

shown  in  Figs.  4-41  to  4-44  for  rf  powers  of  25  watts  and  196  watts. 

It  is  found  that  the  region  of  maximum  rate  of  change  of  phase,  i.e., 
maximum  k^,  the  dip  in  the  rf  amplitude  profile,  and  maximum  electron 
current  density  always  occur  at  the  same  radial  position.  Also  corrobo¬ 
rating  these  results  is  the  fact  that  as  ECRH  power  is  reduced,  the 
radial  position  of  maximum  k^_  moves  radially  inward  as  indicated  in 
Figs.  4-45  and  4-46  to  the  same  radial  position  as  the  maximum  electron 

flux  as  shown  in  Figs.  4-26  and  4-27.  These  rapid  radial  phase  changes 

23 

have  been  observed  by  Fisher  and  Gould  in  the  study  of  the  formation 
of  resonance  cones  in  a  plasma.  It  is  interesting  to  note  from  Figs. 
4-39  to  4-44  that  a  k^Q  <  0  excitation  generates  fields  with  a  positive 
kz  in  the  region  of  electron  flux  which  indicates  electrons  are  moving 
in  the  direction  of  wave  motion.  This  is  probably  due  to  the  f3ct  that 
for  k  <  0  the  source  is  exciting  rf  fields  toward  a  region  of  in¬ 
creasing  electron  density  along  z  and  the  wave  will  experience  a  re- 

22  iO 

flection  as  discussed  by  Eriggs  and  Parker  and  Ott,  et  al. 
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Chapter  5 


THEORETICAL  FORMULATION  AND  INTERPRETATION 
OF  EXPERIMENTAL  RESULTS 

5.1  Introduction 

In  this  chapter  we  consider  the  basic  theoretical  formulation 
of  the  wave-plasma  interaction  and  discuss  the  applicability  to  our 
experimental  results.  The  theoretical  section  of  this  chapter  can 
be  basically  divided  into  two  parts;  the  first  deals  with  the  excita¬ 
tion  of  electrostatic  waves  in  an  inhomogeneous,  magnetized  plasma  by 
a  slow  wave  structure,  and  the  second  part  deals  with  the  interaction 
of  the  electrostatic  waves  with  the  plasma.  While  no  attempt  has  been 
made  to  give  a  complete  rigorous  theory  of  these  problems  a  summary 
of  the  basic  theory  and  relevant  results  from  the  literature  will  be 
discussed . 


5.2  Electrostatic  Waves  in  a  Magnetized  Plasma 


The  excitation  of  electrostatic  waves  in  plasmas  by  using  a  slow 

.  .  ...  u  8,10,12,13 

wave  structure  has  been  studied  extensively  in  the  literature 

An  axial  electric  field,  is  imposed  near  the  plasma-vacuum  interface 

using  the  exciting  structure.  The  ph-'ing  of  the  structure  is  such  that 

•w/kzo<<c  and  hence  an  electrostatic  approximation,  E*-~J  is  assumed. 

The  excitation  of  waves  in  the  frequency  ranpe  .2.  <<  j  <<  d  <<  2  is  dis- 

.  _  '  l  pe  e 

cussed  and  a  time  dependence  e  is  assumed. 

31 

Following  the  notation  used  by  Sti;:.  ,  Maxwell's  equations  can  be 
cast  into  the  following  form. 


V  x  E  ■  iwu  H  (5.1) 

-  o* 

V  x  H  »-iioe  K  •  E  (5.2) 

o  = 

where  K  is  the  cold-magnetized  plasma  dielectric  tensor  defined  by 
Eq.  (5.3) 


Assuming  only  x  and  z  variations  and  taking  the  divergence  of  Eq 
(5.2),  $  can  be  determined  from 


f-  K  +  f  K  !*•  -  0 

dX  XX  oX  dz  ZZ  3Z 


In  the  frequency  range  of  interest,  i.e., 
K  are  approximated  by 


0)  . 

K  «s  1  -  -2~ 

xx  2 

ui 


and 


«  1  — 2®. 


n  <<  u>  <<  w  <<  ft  k 
i  pe  e  xx 


(5.8) 


(5.9) 


I 

I 
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It  should  be  noted  that  K  and  K  are  both  functions  of  x  and  z 

XX  zz 

since  the  electron  density  is  inhomogeneous  in  x  and  z.  Fig.  5-1 
shows  the  geometry  of  our  model.  The  model  in  Fig.  5-1  is  not  strict¬ 
ly  valid  for  our  experimental  results  since  the  model  assumes  a  con¬ 
stant,  uniformly  aligned  density  gradient  at  an  angle  y  with  respect 
to  the  x-axis  and  experimental  results  indicate  a  Gaussian  radial  den¬ 
sity  profile.  However,  since  the  wave  is  found  not  to  propagate  deep 
into  the  plasma  we  will  assume  a  local  constant,  uniformly  aligned 
density  gradient  near  the  periphery  of  the  plasma  column.  Fig.  5-1 
also  defines  a  new  coordinate  system,  (p,n).  The  n  and  o  axes  are 

perpendicular  and  parallel  to  the  density  gradient  respectively  and 

e 

0  ■  constant  defines  surfaces  of  constant  density.  The  new  coordinates 
are  defined  by 


x  «  n sin y +  ocos y  (5.10) 

z  =  n  cos  y  -  psinY  (5.11) 

and  equivalently 


0  z 


Fig.  5-1  Geometry  of  the  Model 


P  *>  x  cos  y  -  z  sin  y 
n  “  xsiny  +  z  cos  y 


(5.12) 

(5.13) 


Eq.  (5.7)  must  in  general  be  solved  in  conjunction  with  the  boundary 
condition  defined  by  the  potential  field  associated  with  the  slow  wave 
structure  as  in  Eq.  (5.14) 


$(x  *  0,2)  -  4>  e 

Q 


ik  2 
20 


(5.14) 


where  is  the  peak  of  the  rf  potential  applied  between  adjacent 
coupler  plates.  Eq.  (5.7)  can  be  rewritten  in  the  (p,n)  coordinate 
system  bv  using  the  chain  rule  and  using  Eas.  (5.12)  and  (5.13)  to 
transform  the  x  and  z  derivatives.  The  advantage  in  doing  this  is  that 
in  the  (o,n)  coordinate  system  the  components’ of  the  dielectric  tensor 
are  only  functions  of  p. 


3$  m  _3p  j)J>  _3_o 
3x  3p  3x  3n  3x 


3$  ,  .  3i> 

cosy-^+  siny^ 


(5.15) 


102 


103 


Dividing  by  cos  y  and  grouping  terras  yields 

[(K  +  K  tan2  y)  ^-1  +  ~  f(K  +K  tan2  y)  -jp-1 

3p  l  xx  zz  3p  J  3n  [  zz  xx  anj 

+  x~  f(K  tan  y  -  K  tan  y)  +  —■  f(K  tan  y  -  K  tan  y)  ] 

2p  L  xx  zz  3nJ  3n  L  xx  zz  3p  J 


(5.18) 

0 


Near  x»0,  k  » k  and  k  :0  for  small  y.  Since  at  x*0  a  wave  is 
n  zo  o  ik  z  ik 

launched  with  the  dependence  e  2°Z  ;  e  ""  and  the  ^electric  tensor 
is  invariant  along  n  we  assume  plane  wave  solutions  of  the  form 


ik  n 

$(p,n)  -  C(a)e  n 

and  ~  ~ ^  ik  . 

3n  n 


(5.19) 


Using  the  plane  wave  solution  indicated  in  Eq .  (5.19)  and  noting  that 

K  and  K  are  independent  of  n  Eq.  (5.18)  can  be  rewritten  as 
XX  z  z 


—  f(K  +  K  tan2y)|£  ]  -  k2  [k  +K  tan2y]  C  + 

3p  [  xx  zz  '3p  J  n  l  zz  xx  'J 

ikn  I^T  ~K  Hanyl  ?+  f(K  -K  )tanyl-p-i-0 

n  f  3p  |_  xx  zz  J  [_  xx  zz  J  3o  j 
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Following  the  notation  used  by  Briggs  and  Parker  we  define 


(5.20) 


(5.21A) 


pn 


(K  -K  )tanv 
zz  xx 


K  »  K  +K  tan  y 
T)  zz  xx 


and  Eq.  (5.20)  is  written  as 


(5.21B) 

(5.21C) 
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K  -  ik 
3p  o  3p  n 


(t~  K  +  K  £  -  k2K 

l 3p  pn  pn  3p  J  nr 


(5.22) 


In  weakly  inhomogeneous  systems,  that  is,  for  large  k  (small  wave- 

P 

length  in  o-direction  compared  with  the  scale  length  in  p-direction) , 
it  is  possible  to  write  a  local  spatially  dependent  dispersion  relation 

assuming  5  -  exp  £ij  kp(p')dp'  j.  Noting  that  ■^-•*-ikp,  we  have  from 

Eq.  (5.22), 


K  k2  +  2k  K  k  -  k2K  -  0 

p  p  n  pn  p  n  n 


(5.23) 


Since  k^  is  fixed  from  the  boundary  condition  at  x-0,  k^  can  be  deter¬ 
mined  from  Eq.  (5.23)  as, 


— k.  rv  K.  .  , 

k  -  -■£•-££  +  (K2  -  K  K  )1/2 

p  K  -  K  pn  p  n 

p  p 

Using  Eqs .  (5.21),  Eq.  (5.24)  can  be  written  as 


(5.24) 


kp  ’  f  KKon  -  (-KXXK„)1''2]  (5.25) 


If  y«0,  i.e.,  Vn  JB  ,  k  ■  k,  k  *  k  ,  K  ■  K  ,  and  K 
e-^-o  p  x’  n  z  p  xx’  < 


0,  yield- 


te)1 


(5.26) 


where  the  plus  sign  was  chosen  to  be  consistent  with  inwardly  plane 
waves . 

Eq.  (5.26)  i3  the  usual  result  obtained  for  a  density  gradient  per¬ 
pendicular  to  the  background  magnetic  field.  From  Eq.  (5.25)  it  is  seen 

that  there  are  propagating  waves  when  K  and  K  are  opposite  in  sign. 

XX  z  z 

This  corresponds  to  the  condition  that  w  ,  <  ui  <  u>  .  Hence  waves  can  be 

pi  pe 

launched  from  the  low  density  periphery  of  the  plasma  column  into  re¬ 
gions  of  increasing  density,  k  is  determined  from 

z 


» 
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k  *k  cos  y  -k  sin  Y  (5.27) 

Z  T\  (5 


Using  k  defined  from  Eq.  (5.25)  and  the  fact  that  k  *k  ,  k  can 
p  n  zo  z 

be  written  as 


k  cos 
zo 


Y-tt-  sin  y  f-k  K  +  |k  | (1 + tan2?) (-K  K  )1/2  1  (5. 

K  L  zo  pn  -  zo 1  xx  zz 


28) 


Similarly 


and 


k  *  k  sin  y  +  k  cos  Y 
x  n  o 


(5.29) 


k 

x 


k  sinY+n- 

ZO  K 

p 


COS  Y 


[- 


k  K  + 
zo  pn  - 


zo 


(l  +  tan2Y)(-K  K  )1/2| 
xx  zz  j 


(5.30) 


Eq.  (5.24),  or  equivalently  Eqs.  (5.28)  and  (5.30)  indicate  two  possible 

mode  structures.  These  modes  and  their  properties  have  been  discussed 

22 

in  detail  for  lower  hybrid  waves  by  Briggs  and  Parker  and  by  Maggs  and 
32 

Morales  for  electrostatic  modes  in  a  magnetized  plasma  with  a  longi¬ 
tudinal  density  gradient.  We  present  here  a  brief  summary  of  the  prop¬ 
erties  of  these  modes.  From  Eq.  (5.24)  we  define  these  modes  as 


-k  K  |k  |  _  .  ,, 

k  -  — ~  -P-  +  -jA-  (K2  -  K  K  )1/2  (5.31) 

d+  K  K  pn  P  n 

p  p 

-k  K  |k  j  .  , 

k  -  -  I!  °-r|'  -  (K2  -  K  K  )1/2  (5.32) 

p-  k.q  pn  o  n 

The  resonance  layer  is  located  at  the  singular  point  of  the 

differential  equation  defining  £  in  Eq.  (5.22).  The  singular  point 

22 

is  determined  from  K  -*■  0 


K 

o 


K  + 
xx 


K  tan2y 
zz 


0 


(5.33) 
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The  density  at  the  resonance  layer  can  be  determined  from  Eq .  (5.33) 
as 


Solving  for  the  critical  density  yields 


Vo  2  /  1  +  tan^Y 
“c-  2  “  (  m  , 

\1+'5_tanY> 


(5.34) 


From  the  experimental  density  profile  given  in  Eq.  (4.2) 


2  .2 

-r  /R  2. 

n  -  n  e  (1  -7-) 

O  L 


(4.2) 


Evaluating  jj  and  -^,the  angle  the  density  gradient  makes  with  the  r- 
axis  is  given  by 


tan  y  ■ 


R 


2 


2(L  -  z)r 


(5.35) 


From  Eq.  (5.35)  it  is  apparent  that  the  angle  y  is  a  function  of  r  and 

z  which  severely  complicates  the  analysis.  If  however,  we  assume  that 

the  critical  density  occurs  at  the  region  of  maximum  current  drive 

(r-3.8  cm)  and  evaluate  Eq.  (5.35)  at  z»-3.81  cm  it  is  found  that  y 

can  be  locally  approximated  by  y*1-58*.  Inserting  experimental  param- 

8  “3 

eters  into  Eq.  (5.34)  yields  n£  :  1.7x10  cm  comparing  favorably  with 

typical  values  of  electron  density  at  these  coordinates  ranging  from 

8  -3  8  -3 

approximately  3.5x10  cm  to  6  x  10  cm  as  determined  from  Langmuir 
probe  measurements  (see  Fig.  4-17).  These  results  are  surprisingly 
accurate  considering  the  approximations  made  and  the  tolerances  in  ex¬ 
perimental  measurements . 
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We  now  consider  a  brief  physical  description  of  the  wave  fields 

in  the  plasma  for  the  two  modes  described  in  Eqs.  (5.31)  and  (5.32). 

First  we  note  that  for  densities  large  enough  such  that  (-K  K  )  <0, 

xx  zz 

i.e.,  w  <  the  wave  fields  are  decaying  (see  Eq.  (5.25))  for  both 

modes.  Noting  that  K<0,  i.e.,w.<u<u  ,  for  k  >0  and  k  <0  and 

p  pi  pe  n  n 

kp  > 0  both  modes  can  propagate.  If  k^  > 0  the  wave  propagates  in  mode 
kp+  and  experiences  a  resonance  as  K  -*0.  If  k^ < 0  the  wave  propagates 
in  mode  kQ  through  the  resonance  layer  and  reflects  from  the  lower  hy¬ 
brid  layer.  It  now  propagates  in  mode  k  and  experiences  a  resonance 
as  K  -*•  0.  If  the  wave  is  launched  into  a  region  of  higher  density 
along  z,  the  wave  number  in  the  z-direction,  k^  as  given  by  Eq.  (5.28) 
can  change  sign  as  the  wave  reflects  from  the  lower  hybrid  layer  and 
approaches  the  resonant-  layer  from  the  high  density  side.  In  Eq.  (5.28), 
near  the  resonance  layer,  and  is  negative;  for  k^  <  0  or  equivalent¬ 

ly  *c20<^>  the  resonance  layer  is  approached  via  mode  kp+,  and  choosing 
the  upper  sign  in  Eq.  (5.28)  we  see  that  the  second  term  has  a  large 

positive  value  and  hence  the  wave  propagates  in  the  +z-direction  even 

— ilc  2 

though  the  source  has  an  e  20  excitation. 

Experimentally  it  is  found  that  large  phase  changes  occur  as  one 
crosses  the  resonance  layer.  As  the  rf  probe  is  moved  toward  more  posi¬ 
tive  z,  i.e.,  lower  electron  density, it  is  found  that  the  phase  change 
is  not  as  large,  and  in  fact  for  some  radial  position,  k^  does  change 
sign  as  indicated  in  Figs. 4-40,  4-42,  and  4-44.  Also,  as  the  electron 
density  is  reduced  by  reducing  ECRH  power  the  region  of  large  k^  moves 
radially  inward  to  a  new  location  where  the  electron  density  is  the  same 
as  previously.  More  about  the  interpretation  of  physical  results  will 
be  discussed  in  section  5.3. 

5,3  Interaction  of  Electrostatic  Waves  with  Plasma 

Even  though  a  large  amount  of  experimental  data  was  compiled  on 
the  electron  dynamics  in  the  presence  of  the  externally  excited  slow 
electrostatic  wave,  it  is  difficult  to  interpret  the  experimental  re¬ 
sults  because  of  the  many  competing  effects  in  the  plasma.  Hence  we 
will  briefly  describe  several  possible  physical  interpretations  of  the 


* 
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experimental  results.  We  also  summarize  the  pertinent  results  from 
the  literature. 

First  we  consider  the  modification  of  the  electron  distribution 
function  due  to  an  electrostatic  wave  in  the  plasma.  Starting  with 
the  collisionless  Vlasov  equation  with  an  electrostatic  perturbation 
♦  (r.t)  -ei(^-wt) 


[—  +  v  •  Vj-  —  E  •  V  f)-  0  (5.36) 

'3t  -  m  ~  v  ' 

e 

where  E«-7$,  it  is  possible  using  quasilinear  theory  to  recast  Eq. 

(5.36)  into  the  quasilinear  diffusion  equation  by  using  a  second  order 

33 

perturbation  technique 


~  m  7  .  d  •  V  f 

3t  v  =v  v 


(5.37) 


The  velocity  space  diffusion  tensor,  is  defined  by 


r  W  /  kk  \ 

?v *  -sf-  \  iTr-7-w)("^J 


dk 


(5.38) 


where  u  ■  <n  +iu, 
r  i 


and  Wk  is  the  spectral  energy  density  of  the  electrostatic  field.  The 
average  electric  energy  density  of  the  electrostatic  field  is  given  by 


<U_>  -  i  e <  1 E 1 2 >  -  [  W.  dk 
E  2  11  j  k  * 


(5.39) 


and 


W  m 

k 


2(2it)  3V 


(5.40) 
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The  Fourier  transform  of  the  electric  field  is  defined  in  Che  usual 
sense  as 


?k(t) 


(2ir)  3 


E(r, t)e 


ik 


(5.41) 


The  spectral  energy  evolves  in  time  according  to 

Wfc(t)  -  2Wi(k,t)Wk(t)  (5.42) 

We  are  primarily  concerned  with  the  distribution  function  in  axial 
velocity  and  hence  we  only  consider  the  one-dimensional  form  of  Eq. 

(5.37) 


2  W 

-c  2io  .  f  k 

11  =  _E£  _JL _ 1 _  dk  JLL 

3t  n  3v  i(k  v  -  w)  z  3v 
O  Z'  ZZ  z 


(5.43) 


The  axial  electron  kinetic  energy  density  is  given  as 


..12 

LKE  ‘  2  noVz 


(5.44) 


Multiplying  the  one-dimensional  diffusion  equation  by  the  axial  electron 
kinetic  energy  density  and  integrating  over  axial  velocities  yields  a 
conservation  of  energy  relation. 


3  /  1 

3t  '  2  nome 


v^f (v  ,t)dv 
z  z  z 


)  -  ~  fw 

'  3c  k 


(t)dk 


(5.45) 


Eq.  (5.45)  indicates  that  as  the  electrostatic  wave  energy  decreases 
by  some  damning  mechanism  the  electron  kinetic  energy  must  increase  by 
a  corresponding  amount. 

The  exciting  slow  wave  structure  used  in  our  experiment  excites 
a  relatively  broad  spectrum  in  of  electrostatic  waves  since  the  ex¬ 
citing  structure  is  only  one  wavelength  long,  i.e.,  Lg"2n/k2o>  where 
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Lg  is  the  axial  extent  of  the  slow  wave  structure.  The  spectral 
energy  density  associated  with  the  exciting  structure  can  be  calc¬ 
ulated  using  Eqs.  (5.40)  and  (5.41) 

,Ls  i(k  -k  )z 

£  e  20  2  dz  (5.46) 
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Evaluating  yields 
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<k  -k  r 

zo  z 


sin2  \  (k  -k  )L 
2  zo  z  s 


(5.47) 


(5.48) 


It  is  easily  seen  that  the  main  peak  of  the  spectral  energy  is  con¬ 
tained  in  the  wave  number  range  0<kz<2kzo  which  corresponds  to 
electron  energies  greater  than  approximately  0.67  eV  for  kzo *  .21  cm  *. 
Hence  the  bulk  of  the  electron  energy  distribution  should  be  modified 
by  the  presence  of  the  electrostatic  wave.  This  results  in  an  increase 
in  the  average  axial  electron  kinetic  energy  which  manifests  itself  as 
an  apparent  increase  in  the  electron  temperature  (see  Figs.  4-32  and 

4.35).  This  energy  is  an  energy  associated  with  the  oscillations  of 
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the  electrons,  i.e.,  it  is  not  thermalized.  It  has  been  shown 
that  the  effective  electron  temperature  associated  with  the  plasma 
electrons  is  given  by 


T 


eff 


(5.49) 
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where  T^^  and  Tg0  are  measured  in  eV.  Experimentally  it  is  found 

that  the  effective  temperature  increases  strongly  only  near  the  outer 

regions  of  the  plasma,  i.e.,  near  the  region  of  maximum  electron  flux 

and  radial  phase  change.  These  results  are  consistent  with  the  results 
13 

of  Decyk,  et.al.  in  their  observation  of  electron  surface  heating  us¬ 
ing  a  phased  array  of  capacitor  plates  to  couple  electrostatic  energy 
into  the  plasma. 

We  will  now  discuss  a  simple  physical  interpretation  of  the  rf- 
induced  axial  electron  flux  near  the  plasma  surface  based  on  the  pon- 
deromotive  force  on  electrons  in  an  inhomogeneous  electrostatic  poten¬ 
tial  in  a  background  magnetic  field  Bq. 

3v  _P 

—  +  v  •  7v  ■  —  (E  +  vxB  )  (5 .50) 

3 1  -  -  m  -  -  -o 

e 

and 

E»-7(j>  (5.51) 


The  oscillating  electric  field  has  components  in  both  the  x  and  z 
directions.  Noting  chat 

v  •  Vv  -  |  7(v2)  -  v  x  7  x  v  (5.52) 

the  z-component  of  Eq,  (5.50)  can  be  written  as 
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(5.53) 


It 

3y 


is  assumed  that  there  is  no  variation 

-  0. 


in  the  y-direction. 


i.e. 
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From  Eq.  (5.53)  it  is  desired  to  determine  the  second  order  time 
averaged  nonlinear  force  on  the  electrons  along  the  background  mag¬ 
netic  field.  To  do  this  we  compute  the  first  order  velocities  assoc¬ 
iated  with  the  electrostatic  wave  using  the  linearized  equation  of 
motion  and  then  computing  the  z-component  of  the  ponderomotive  force 
by  time  averaging  Eq.  (5.53).  Writing  the  linearized  equations,  i.e. 
neglecting  the  convective  term  v  -  7v  from  Eq.  (5.50)  and  assuming  a 
time  dependence  e  lu)t 
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(5.54) 

(5.55) 

(5.56) 


Solving  for  v^  yields  from  Eqs .  (5.55)  and  (5.56) 
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Time  averaging  Eq.  (5.59)  yields 
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Noting  that 
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3  2,  _ x 

—  v  *  7  v  — — 
3z  x  x  3z 


and 


,  ,  3v 

3  2,  v 

—  v  *  2v  -T-— 
3z  y  v  3z 


Eq.  (5.60)  can  be  written  as 


3v  .  ,  ,  3v 

_ z  >  _  .  _  1  3  _V  _ 5.  > 

3t  2  3z  Vz  x  3x 


Using  the  results  from  Eqs.  (5.54)  to  (5.57)  we  have 
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The  ponderomotive 
as 


force  density  on  the  electrons. 
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F  can  be  written 
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(5.68) 


Hence , 
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This  result  has  been  obtained  bv  Morales  and  Lee  in  the  study  of 
the  nonlinear  f ilamentation  of  lower-hybrid  cones.  Since  the  elec¬ 
trons  are  stronglv  magnetized  the  first  term  in  Eq .  (5.69)  dominates. 
This  vields 
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e  u) 

o  pe 


(5.70) 


The  wave  field  possesses  strong  spatial  variations  in  the  axial  elec¬ 
tric  field  and  the  electrons  experience  a  ponderomotive  force  given  by 
Eq.  (5.70)  and  the  electrons  are  forced  out  of  the  region  of  large  rf 
potential.  Experimentally  this  is  seen  as  an  increase  in  the  axial 
electron  flux.  The  electrons  leave  the  high  field  region  leaving  a 
net  positive  charge  produced  by  the  ions  which  produces  a  large  elec¬ 
tric  field.  The  ions  accelerate  out  oF  this  region  to  maintain  auasi- 

neutralitv.  This  leads  to  a  localized  density  depression  in  the  high 

35  36  37 

field  region  as  discussed  bv  Morales  and  LeeJ  ,  Oekelman  and  Stenzel  ’  ' , 

38  1 3 

Wilson  and  Wong  ,  and  Oecvk,  et.al. 

5.4  Conclusions  and  Proposal  for  Future  Work 

In  this  work,  we  have  presented  experimental  results  on  the  inter¬ 
action  of  electrostatic  waves  excited  by  a  slow  wave  structure  with  the 
plasma.  It  is  found  that  large  radial  phase  changes  exist  near  the  lower- 
hybrid  resonance  layer  and  the  wave  launched  from  the  slow  wave  structure 
can  experience  a  reflection  when  the  wave  is  launched  along  z  into  the 
high  density  region  of  the  plasma.  Large  increases  in  the  effective 
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axial  electron  temperature  are  found  near  the  plasma  surface  consistent 
with  quasilinear  theory.  Increases  in  the  axial  electron  flux  can  be 
attributed  to  the  large  ponderomotive  force  on  the  electrons  which  can 
lead  to  ion  acceleration  and  a  density  depression  due  to  the  transient 
space  charge  field  as  the  electrons  leave  the  high  rf  field  region. 

More  detailed  phase  measurements  should  be  made  in  the  future  to 
obtain  a  greater  understanding  of  the  wave  structure.  This  should  in¬ 
clude  azimuthal  phase  measurements  so  that  a  model  may  be  developed 
in  cyclindrical  geometry.  Also  the  ion  dynamics  should  be  studied  by 
using  the  electrostatic  energy  analyzer  to  measure  the  ion  distribu¬ 
tion  function  to  verify  the  self-consistent  acceleration  of  ions  due 
to  the  ponderomotive  force  on  the  electrons.  In  addition  the  mechanism 
associated  with  plateau  formation  in  the  distribution  function  as  well 

as  the  saturation  of  the  rf-incuded  electron  flux  for  a  k  <0  excita- 

zo 

tion  should  be  investigated. 
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